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1 Introduction 

After Elie Cartan, who in 1922 introduced differential geometric methods into the calculus of 
variations using the celebrated Cartan form to formulate the variational principle [10], Theophile 
Lepage in his pioneer work [73, 74] initiated the study of the calculus of variations for field theo¬ 
ries (multiple integrals) as a theory of differential forms and their exterior differential modulo 
contact forms (congruences). The main advantage of such an approach evidently is a possi¬ 
bility of framing the calculus of variations within the de Rham sequence of differential forms 
on a manifold and thus using the de Rham cohomology as an effective tool to study relations 
between the image and kernel of various differential operators appearing in the calculus of 
variations. 

The attempt to understand global structures of the calculus variations motivated the use of C- 
spectral sequences, and gave rise to the constructions of variational sequences and bicomplexes. 
This direction was initiated in 1950’s by Dedecker [14, 15, 16, 17, 18], and continued in the work of 
Horndeski [41], Gelfand and Dikii [31], Olver and Shakiban [79], Dedecker and Tulczyjew [19], 
Tulczyjew [96, 97, 98], Takens [93], Anderson [2], Anderson and Duchamp [3], Dickey [20], 
Vinogradov [99, 100, 101], and others. A main motivation for this research was 

• to find solution of the inverse problem of the calculus of variations which means to find 
local and global constraints for a system of ordinary or partial differential equations to 

*This paper is a contribution to the Special Issue on Analytical Mechanics and Differential Geometry in honour 
of Sergio Benenti. The full collection is available at http://www.emis.de/journals/SIGMA/Benenti.html 
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come from a variational principle as Euler-Lagrange equations, and to construct a Lag- 
rangian, 

• to determine the local and global structure of variationally trivial Lagrangians, i.e., Lag- 
rangians giving rise to identically zero Euler-Lagrange expressions (the kernel of the Euler- 
Lagrange operator), 

• to understand the form and origin of different “canonically looking” expressions appearing 
in the calculus of variations, as for example 


dL 

dy^ 


, dL ^ ^ 

^3 QyU + ^3^k 


dL 


dEa dE, 

H ^ \j) dql ’ 

d K—J 


dL r dL 


dL 


■ip 

s-j-l 

E 

1=0 


ijp 


dy 


jkp 




j + A S dE, 
) dt^ dq 


\ I 


U 

'j+i+i 


• to understand local and global phenomena typical for variational problems given by lo¬ 
cal data, as, for instance examples of global Euler-Lagrange forms which do not ad¬ 
mit a global Lagrangian; the best known examples are the Galilean relativistic mecha¬ 
nics, or Chern-Simons field theories (where a global Lagrangian exists but is not gauge- 
invariant). 

The theory was first formulated on infinite order jets. After several attempts to obtain an 
alternative construction using finite order structures (see, e.g., [3]), the finite order variational 
sequence was discovered by Krupka 25 years ago [49]. The inspiration originated from the 
Lepage’s idea of a ‘congruence’, indicating that there should be a close relationship between the 
Euler-Lagrange operator and the exterior derivative of differential forms, and resulted in the 
concept of a finite order (cohomological) exact sequence of forms, in which the Euler-Lagrange 
mapping is included as one arrow. 

Compared with the other approaches, the concept of Krupka’s variational sequence is con¬ 
ceptually very simple due to the following reasons: 

• it uses only finite order jets, 

• it appears as a quotient sequence of the de Rham sequence. 

This makes it quite easily accessible to potential users, since factorization is a natural, 
straightforward, commonly understood, and well-mastered technique. Moreover, the factori¬ 
zation procedure can be incorporated into intrinsic operators which allow to avoid tedious and 
technically difficult calculations. 

It should be stressed, however, that here is an additional benefit of the setting on finite order 
jets: Namely, it provides precise results on the order of the local and global objects of interest. 

The aim of this paper is to provide the reader a comprehensive source of the current theory 
of the finite order variational sequences as it stands 25 years after its beginnings. While for the 
variational bicomplex theory there is a standard reference the (unpublished) book by Ander¬ 
son [2], the results on the finite order variational sequence up to now have not been collected, 
revised and completed in a self-contained article. However, our aim is not merely to provide 
a revised exposition of known results. Apart from putting known results into a unified treat¬ 
ment, we add a number of original results in order to develop the theory and make it complete. 
Our aim is also to stimulate the use and applications of the variational sequence in physics. As 
a motivation we mention a few applications, some of which are interesting and not widely known 
or even are original here. 

The paper is organized as follows: Section 2 introduces the variational sequence as a quotient 
sequence of the de Rham sequence, developed by Krupka [49, 51, 52, 53] (see also [36, 37, 54, 55]). 
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Section 3 is the core of the paper, and deals with the problem of representation of the va¬ 
riational sequence by differential forms. In the variational sequence, the objects are classes of 
local differential forms. However (as pointed out already by Krupka) the classes can be repre¬ 
sented by (global) differential forms: Lagrangians or Lepage n-forms (where n is the number 
of independent variables in the variational problem), Euler-Lagrange forms, Helmholtz forms. 
This property initiated the search of representation operators and representations with help of 
differential forms. We present two different representation sequences: the Takens representa¬ 
tion sequence based on the so-called interior Euler operator [2] which assigns to every class in 
the g-th column of the variational sequence a source form of degree q (a Lagrangian if g = n, 
a dynamical form if (7 = n-|-l, a Helmholtz-like form if g = n-|-2), and the Lepage representation 
sequence. The latter is a sequence of Lepage forms (Lepage equivalents of source forms), such 
that the morphisms are just the exterior derivatives. The existence of such a sequence clarifies 
the relationship of the variational morphisms to the exterior derivative operator. Moreover, 
it transfers the questions on the structure of the kernel and image of all the variational mor¬ 
phisms just to application of the Poincare lemma. A complete representation of the variational 
sequence by source forms (of any degree) was given by Musilova and Krbek [43, 44, 45] and 
independently by Krupka and his collaborators [61, 63, 103]. Next, in Section 3 we recall the 
concept of Lepage form ([46, 61, 64, 72] and introduce the Lepage representation of the varia¬ 
tional sequence. As a new result concerning Lepage n-forms, we present an intrinsic formula 
for the Cartan form (Theorem 3.8, formula (3.6)). Remarkably, our formula explicitly shows 
that (and explains why) the higher-order Cartan form generically is not global (this fact was 
a surprise when discovered and demonstrated by a counterexample in 1980’s [40], and up to 
now has not been well understood). In the subsection on Lepage forms of higher degrees (and 
any order) we present for the first time properties of general Lepage forms, and namely we 
give a general intrinsic formula for a Lepage (n -|- A:)-form of any order. We also find a higher 
degree generalization of the famous Cartan form. The techniques based on the use of the in¬ 
terior Euler operator and the residual operator not only lead to elegant formulas and easy 
proofs; they also enable to avoid lengthy, tedious, and sometimes even impossible calculations 
in coordinates. The results then help us to find a complete correct Lepage representation of 
the variational sequence. The last remark in Section 3 concerns the variational order of ob¬ 
jects in the variational sequence. We explain how due to the finite order setting it is possible 
to make (local and global) conclusions on a precise order of Lagrangians and other objects 
in the variational sequence, distinguishing it from alternative constructions on infinite order 
jets. 

In Section 4 we are concerned with the Lie derivative in the variational sequence. We explore 
the important property that the Lie derivative of differential forms with respect to prolongations 
of projectable vector fields preserves the contact structure. This suggests that a Lie derivative 
of classes of forms, a variational Lie derivative, can be correctly defined as the equivalence class 
of the standard Lie derivative of forms and represented by forms. New results in this section are 
then contained in subsequent theorems generalizing the first variation formula to forms of any 
degree. 

The last section of the paper is devoted to selected applications of the variational sequence 
theory. We recall the fundamental well-known applications to the inverse variational problem 
and to variationally trivial Lagrangians. Apart from that we emphasize some other possible 
applications which have not yet been fully explored - like, for instance, properties of Helmholtz 
forms. We include also examples of variational problems defined by local data, and on sym¬ 
metries and conservation laws, emphasizing cohomological questions and local and global phe¬ 
nomena. 

Eor a reader interested in a wider context of the theory of variational sequences, we refer to 
the survey works [54, 102], and to the book [55]. 
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2 The variational sequence 

2.1 Jet bundles and contact forms 

Fibred manifolds and their jet prolongations represent a convenient mathematical framework 
for mechanics and field theories, suitable for investigation of Lagrangian systems of different 
orders, and one or many independent variables within a sufficiently general and unified geometric 
framework. For an excellent exposition of the jet bundle geometry we refer to the book [91]. 

In what follows, we shall consider a smooth fibred manifold vr: V —>■ X, with dim = n and 
dim V = re + m, and its r-jet prolongations vr^: J^Y —)• X, where r > 1. For re = 1 the mani¬ 
fold V is a space of events for mechanical systems of m degrees of freedom, and local sections 
of vr are graphs of curves, so usually X = M (meaning time, or just a parameter for the curves). 
If re > 1 then X has the meaning of a physical space or spacetime, and local sections of tt are 
physical fields over the manifold X. Then J^Y is a manifold of points jjjy, the equivalence 
classes of local sections 7 of vr with the same value at x and the same partial derivatives at x 
up to the order r. Thus J’^Y serves as a domain for functions, differential forms, vector fields, 
and other objects depending on higher derivatives (up to the r-th order). 

Due to the affine bundle structure of vr^+i^r-: J''~^^Y —)• J'’X, we have a natural splitting 
J^Y XjT-iy T*J'~~^Y = J^Y Xjr-iy {T*X © V* J'^~^Y) which induces natural splittings into 
horizontal and vertical parts of projectable vector fields, forms, and of the exterior differential 
on rY. 

A differential g-form oj on J^Y is called contact if = 0 for every local section 7 

of vr. Contact forms on J^Y form an ideal in the exterior algebra, called the contact ideal of 
order r. The contact ideal is generated by contact I-forms and their exterior derivatives [50]. 
A differential q-form to on J'^'Y is called horizontal, or 0-contact, if = 0 for every vector 
field ^ on J’’X, vertical with respect to the projection onto X. 

From the definitions it follows that every g-form for > re is contact, and that every form uj 
on rY, if lifted to r+^Y, is a sum of a unique horizontal and contact form. This splitting can 
be further refined if the concept of a k-contact form for k > 1 is introduced. A contact form uj 
is called fc-contact if for every vertical ^ the contraction is {k — I)-contact. Then, as proved 
in [48], every g-form a; on J’’X, g > I, admits the unique decomposition 

K+i,t^ = ( 2 - 1 ) 

k=0 

into a sum of /s-contact forms, 0 < k < q. For the horizontal (0-contact) operator often the 
notation h = po is used. The form is called the k-contact component of uj. Note that if 
g > re + I then u is contact, and it is at least (g — n)-contact, i.e., the contact components 
Po,Pi,. ■ ■ ,Pq-n-i of UJ are equal to zero. If, moreover, Pq-nOJ = 0, we speak about a strongly 
contact form [49] . 

An important special case of the decomposition (2.1) concerns closed forms and their local 
primitives, obtained due to the Poincare lemma, by means of the Poincare homotopy opera¬ 
tor V. Since every form uj can be locally expressed as a; = Vduj + dVuj, equation uj = dp has 
a local solution p = Vuo. Applying, however, the decomposition into contact components we 
can see that there arises a modified homotopy operator, denoted by A, and called the contact 
homotopy operator [48], which is adapted to the decomposition (2.1) and defined by Apou = 0, 
ApkUJ = Pk-i'Puo, k > 1. It satisfies ^T*_^_l^^uJ = Ad{^T*_^_l^^uJ) + dA{^T*_^_l^^uJ), and moreover, if 
(locally) UJ = dp, then 

<? 

'^r+l,rP ~ •^{'^r+l,r^) ~ Apk+l^i 

k=Q 
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hence pkP = Apk+iOJ. Compared to R, the operator A concerns vertical curves (curves in the 
fibres over X) only. 

Notice that the decomposition (2.1) applied to uj = dp induces a splitting of the exterior 
derivative d to ^d = dn + dy where the horizontal derivative dn = hd, and the vertical 
derivative dy is the contact part of ^d. 

The splitting of vector fields then arises with help of the total derivative operator. If we 
denote local fibred coordinates on Y, and {x^,yj) where J is a multiindex, 0 < | J| < r, 

the associated coordinates on then the z-th total (formal) derivative takes the form 


di 


JL 


+ X] yji 

|J|=0 


d 

dyj' 


Considered as a vector field along the projection Tir+i^r, the total derivative induces a splitting 
of vector fields to a horizontal and vertical part as follows: Letting 


P = A— + V — 

|J|=0 


be a vector field on J^Y, in order to obtain the splitting ^ + Cv; we put 

d 


in = Cdi, iv = i-iH = ^ (?j- y'jiC) 

|j|=o 


dy° 


Note that both and ^y are vector fields along the projection 'Kr+i,r rather than ‘ordinary’ 
vector fields on J'^Y. If H is a projectable vector field on Y, we write = (J’’H)y (this can 

be viewed as a definition of prolongation of Hy which is a vector field along the projection tti^q)- 
Recall that C{Trr), the contact distribution or Cartan distribution of order r, is locally anni¬ 
hilated by contact 1-forms 


wj = dyj - yjjdx ^, 1 < u < m, 

or, equivalently, locally generated by the vector fields di along TTr,r-i and d/dyj where | J| = r. 
The Cartan distribution on J’^Y is not completely integrable, hence the ideal generated by 
contact 1-forms of order r is not closed. Its closure is then the contact ideal of order r (the ideal 
of all contact forms on J'^Y). 

In what follows, we shall need on J'^Y the ideal generated by lifts of contact 1-forms of 
order 1, whose elements take the form A p where uj G C*^(7ri) (the annihilator of C(7ri)) 

and p is of order r. Since locally they are expressed as A p^, where = dy^ — yjdx^, 
1 < a < m, are the basic local contact 1-forms of order 1, we say that forms belonging to this 
ideal are -generated. 

In the calculus of variations and the theory of differential equations on fibred manifolds, the 
fundamental role is played by the sheaf x horizontal (0-contact) n-forms on J^'Y, the 
elements of which are called Lagrangians of order r, and the sheaf i -y of 1-contact (n -|- 1)- 
forms on horizontal with respect to the projection onto V, the elements of which are called 
dynamical forms of order r. By means of variation of the action defined by a Lagrangian A 
one obtains a distinguished dynamical form Ex, called the Euler-Lagrange form of A [46]; the 
components of Ex in every fibred chart are the Euler-Lagrange expressions of the Lagrangian A. 
Remarkably, if A is of order r then Ex is of order < 2r. 

The mapping assigning to every Lagrangian its Euler-Lagrange form is called the Euler- 
Lagrange mapping. 
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A dynamical form E is called globally variational if there exists a Lagrangian A such that 
(possibly up to a jet projection) E = E\. E is called locally variational if every point in the 
domain of E has a neighbourhood U where E is variational [48]. 

Dynamical forms are often called source forms, as originally introduced by Takens [93]. How¬ 
ever, we shall follow Krupkova and Prince [70] who extended the concept of source form to forms 
of all degrees q = n + k where A: > 1: 

Definition 2.1. Let A: > 1. By a source form of degree n + A: we shall mean a cu'^-generated 
A:-contact (n -|- A:)-form. 

For A: = 1, dynamical forms are w'^-generated 1-contact (n -|- l)-forms, hence, indeed, source 
forms of degree n -|- 1. Source forms of degree n -|- 2 are called Helmholtz-like forms [69]. For 
the calculus of variations the most important examples of source forms of degree n -|- 1 are 
Euler-Lagrange forms, and of source forms of degree n -|- 2 the Helmholtz forms. 

2.2 The variational sequence 

Now we recall the finite order variational sequence as introduced by Krupka in [49]. The concept 
uses the sheaf theory (standard references are, e.g., [ 8 , 104]). However, for a first understanding, 
a concrete use of the variational sequence, and calculations it is not necessary to be familiar with 
all the rather difficult mathematical theory. In place of a ‘sheaf’ of forms in this context one 
can roughly imagine a family of modules of local differential forms defined on the open subsets 
of the corresponding manifold. 

In what follows, we denote by My the constant sheaf over M, and by the sheaf of g-forms 
on HY. We let Aq^^ = {0}, and A^^^ be the sheaf of contact g-forms, if g < n, or the sheaf of 
strongly contact g-forms, if g > n, on HY. We set 

0 ^ = A;, + (dA;_i,,), 

where is the image of by the exterior derivative d, and denotes the sheaf 

generated by the presheaf ciAg_^We note that 0^ = 0 for g > corank C( 7 rr). 

The subsequence 

0 ^ 0 ^ ^ 0 ^ ^ 0 ^ ^ • 

of the de Rham sequence 

0 ^ My ^ AS ^ A^ ^ A^ ^ A^ ^ • 

is an exact sequence of soft sheaves. The quotient sequence 

0 ^ My ^ AS ^ A\/Q'[ AS/0S ^ AS/0S ^ • (2.2) 

is called the variational sequence of order r. The main theorem due to Krupka [49] then states 
that the variational sequence is an acyclie resolution of the constant sheaf My (meaning that 
the sequence is locally exact with the exception of My). Hence, due to the abstract de Rham 
theorem, the cohomology groups of the cochain complex of global sections of the variational 
sequence are identified with the de Rham cohomology groups H^j^Y of the manifold Y. 

It should be stressed that objects in the variational sequence are elements of the quotient 
sheaves A'^/Q'^^, i.e., they are equivalence classes of local r-th order differential q-forms. We 
denote by [p] G the class of p G A^. Morphisms in the variational sequence are then by 

construction quotients of the exterior derivative d. We denote 


8 ,-. AS/0S ^ AS+1/0S+1 
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so that 

£q{[p\) = [dp]^ 

and £q{[p]) = [dp\ = 0 means that there exists [p] G ^q_i/Qq_i such that [p\ = £q-i{\p\) = [dp]. 
If, moreover, = {0} then the class [p] has a global representative. 

The meaning of this abstract construction for the calculus of variations becomes clear in view 
of the following theorem: 

Theorem 2.2 (Krupka [49, 54]). The sheaf is isomorphic with a subsheaf of the sheaf 

of Lagrangians the s/iea /is isomorphic with a subsheaf of the sheaf of source 

forms ^de quotient mapping A])/0[[ —)■ A5)_,_]^/0[)_|_^ is the Euler-Lagrange 

mapping. 

By this theorem, the elements of the quotient sheaf A])/0[[ have the meaning of Lagrangians. 
The kernel of £n then consists of null-Lagrangians (variationally trivial Lagrangians), and its 
image, which, by exactness is the kernel of the next morphism, 

a;;+i/0;+i ^ a;;+2/0;;+2, 

represents variational equations. Therefore £n+i is called Helmholtz mapping. The name refers 
to the famous “Helmholtz conditions” [39], necessary and sufficient conditions for differential 
equations to be variational (as they stand). Remarkably, the Helmholtz morphism, discovered 
within the variational sequence theory, has not been known in the classical calculus of variations. 

The above theorem also suggests the idea that all classes in the variational sequence could be 
representable by differential forms with a clear place in the calculus of variations. So far there 
have been noticed two such representation sequences: the “'source form representation” which we 
call Tokens representation, and the “Lepage form representation”. In the former the sheaves con¬ 
tain familiar differential forms which appear in the calculus of variations, like Lagrangians and 
Euler-Lagrange forms, and their generalizations, like, for example, Helmholtz forms. The mor- 
phisms are then the corresponding variational mappings, like, e.g., the Euler-Lagrange mapping 
sending Lagrangians to Euler-Lagrange forms, or the Helmholtz mapping, sending dynamical 
forms to Helmholtz forms. The latter representation has as objects Lepage forms, which are 
extensions of the “familiar forms” mentioned above by higher contact components in such a way 
that the morphisms are the exterior derivatives, so that in the Lepage form representation the 
variational sequence “morally” (up to the orders of individual subsheaves) becomes a subse¬ 
quence of the de Rham sequence. We shall deal with the representation sequences below. 

In what follows, we shall denote 

Vo^ = AS, V”q = Ay%, 

and the variational sequence of order r will be shortly written in the form 

0 ^ Mr ^ Vr. 

3 Representation sequences 

3.1 The interior Euler operator 

Let us introduce an operator which plays an essential role in the representation theory for 
the variational sequence. It was introduced to the calculus of variations within the variational 
bicomplex theory in [2, 4], and adapted to the finite order situation of the variational sequence 
in [43, 44, 45, 61]. This operator, called interior Euler operator, and denoted by X, reflects in 
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an intrinsic way the procedure of getting a distinguished representative of a class [p] G Aq/ 0 q 
for g > re by applying to p the operator Pq-n and a factorization by 0 ^. 

Let denote local fibred coordinates on Y, and {x'^,yj) the associated coordinates 

on J^Y, where J is a multiindex, 0 < \ J\ < r, J = {ji,... ,jp) with I < ji < j 2 < • • • < jp < n 
and p < r. 

Let k > 1 and q = n + k. Recall that if p G ^n+k^ have 

'^r+l,rP ~ PkP + Pk+lP + ‘ ‘ ‘ + Pk+nP- 
Following [45] we set 

^{p) = A ^ {-l)\-^\dj{d/dy^j\pkp), (3.1) 

|J|=o 

where dj = dj^ ... dj^ for | J| = p. It can be shown by means of the partition of unity arguments 
that this formula defines a global form 2 i(p). 

The operator I: —)■ A^^^ has the following properties: 

(1) for every p G A])_|_^, T{p) is a source form of degree n + k, 

(2) Z{pkp) =Z{p), 

(3) Z{p) belongs to the same class as TT^r+i^rP^ i-®-) '^ 2 r+i,rP ~^{p) ^ 

(4) = Z, up to a canonical projection; precisely, Z^{p) = 2 r+i^{p)^ 

(5) KerX = 0];+,. 

By construction, Z{p) is a A:-contact form, so there is a question about the difference bet¬ 
ween Z{p) and PkP- It turns out that the difference can be expressed intrinsically by means of 
an operator TZ which we shall call residual operator, as follows 

PkP = Z{p) + pkdpkTZ{p) (3.2) 

(up to appropriate canonical projections). By the properties (2) and (4) of Z we can see that 

Z{pkdpkTZ{p)) = 0. (3.3) 

TZ{p) is by definition a local strongly contact {n + k — l)-form. Note that, contrary to Z{p), the 
form TZ{p) need not be unique and need not be global. 

Moreover, Z{p) is a source form, so there is a question about its relationship with other source 
forms in the same class [pj. To clarify this it is useful to introduce the concept of a canonical 
source form as a source form p such that (up to projection) p = Z{p) [70]. 

Proposition 3.1. Every source form is equivalent with a canonical source form {of possibly 
higher order). More precisely, letting p = A po- be a source (re -|- k)-form it holds {up to 
projection) 

p = kZ{p) — {k — l)uj^ AZ{r](j) 


and 

k — 1 

p-Z{p) =pkdpkZi{p) = Apk-idpk-iTZ{r]a)- 

As a consequence, for k >2, 

Z{oj'^ A ri„) = Z{uj'^ A Z{ri„)). 
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Proof. The first assertion is just the property (3) of X. Next two formulas easily follow from 
the definitions of X and TZ and the fact that p is the source form, i.e., rP = PkP- Indeed, 
for p = uj^ /\ rjfj we have 

X{uj^ A = iw'" A ^ A p^)) 

|J|=0 

= -uj‘' A ^{-l)^'^Uj{d/dyj\pty)] = A f\X{p„). 

^ |J|=o ^ 

This formula then yields 

PkdpkTZ{p) = p-X{p) = A {pa -X{pa)) = A Pk-ldpk-lTZ{Pa), 

as desired. 

Finally, applying X to the first formula we get {k — l)X{p) = (k — l)X{uj^ AX{p„)), which 
provides the last formula for k >2. ■ 


For dynamical forms Proposition 3.1 gives the following important uniqueness result [70]. 

Corollary 3.2. For k = 1, every source form is a canonical source form. 

Note that this can be seen also directly from the definition of X: Indeed, if /c = 1 and p is 
a cu'^-generated 1-contact (n+ l)-form then in (3.1) the contractions give horizontal n-forms with 
the same components as those of p, and the summation over | J| reduces to one term, jjj = 0. 
Hence, X{p) = A {d/dy^\p) = p. 

For the description of the image of X we refer also to [63] . 

Example 3.3. We shall show the computations of X{p) and 7^(/o) for a second-order 1-contact 
2-form p in mechanics (n=l,A: = l,r = 2). 

Consider local fibred coordinates (t, q^) on Y and the corresponding associated coordinates 
1 < j < r, 1 < a < m, on J^Y, and denote dt the total derivative operator and 
LjJ = dpj — qj_^_ldt the local basic contact 1-forms, and put = ^o- calculation of the 

image of p by the interior Euler operator is direct, using (3.1). We have 

pip = Adt + A^cof Adt + A dt, 


hence 


r=2 


j=0 


The residual term is then 



(A° - dtAl + dfADoj"^ A dt. 


pidpiTZ{p) = TTg 3P1P - X{p) = {dtAl - djADu)"^ Adt + Alu:f Adt + Aluf A dt. 

Let us compute 'R,{p). Since TZ{p) is a contact 1 -form of order two, it holds 7^(/o) = + 

Clujf + Cluj2 = piTZ{p). The 1-contact part of the exterior derivative is pidTZ{p) = —dtClu^ A 
dt — {dtCl + CDujf A dt — {dtCl + Cl)uj2 A dt — A dt. Comparing this with the formula 

above yields 


cl = o, 


dtC: 


+ Cl = -Al, 


dtCl + Cl = -Ai 


-dtCl = dtAl - dU^ 


From here we immediately see that Cl = —Al and (7° = —Al + dtAl, and that the equation 
-dtCl = dtAl - dlAl is satisfied identically. So finally, 

TZ{p) = {dtAl-Al)u:^-Aluf. 


We can easily verify that vr^gpip = X{p) + pidpiR.{p). 
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3.2 Takens representation of the variational sequence 

Once we have the interior Euler operator, we can obtain a sequence of sheaves of differential 
forms (rather than of classes of differential forms), such that both the objects and the morphisms 
have a straightforward interpretation in the calculus of variations. Such a representation of the 
variational sequence, introduced by Krbek and Musilova (see [43, 44, 45]) basically concerns 
source forms {of all degrees), therefore we call it in honour of Takens the Takens representation. 
The motivation for the representation comes from the following observations: 

• If q < n then for every g-form p of order r, the contact decomposition ( 2 . 1 ) takes the form 

K+l,rP = hp+Plp-\ -hPq/9, 

and 0g = Ag,, + {dAf^_if), where we have sheaves of contact forms. This means that given 
a class [p\ G Vq = Aq/ 0 q, one has pi ~ p 2 iff pi — p 2 = fd + dp where ft is a contact g-form and p 
is a contact {q — l)-form, i.e., 

hpi = hp2. 

Hence, every class [p] is completely determined by a unique horizontal form, and, in particular, 
if g = n, by a Lagrangian. 

Summarizing, for q < n the operator of horizontalization, 
h: Ag —)• Ag(^ C p—!■ hp, 

induces a representation mapping, 

R,: V;^A^+i, Rff[p]) = hp. 

• If q = n + k for k > 1 then for every g-form p of order r, the contact decomposition (2.1) 
takes the form 

T^r+l,rP — PkP + Pk+lP + ' ' ' + PqP, 

and 0q = Ag ,, + {dA'^_i ff, where we have sheaves of strongly contact forms (with the only 
exception of Aj),, which is a sheaf of contact forms). This means that given a class [p] G — 
P^ ~ /^2 iff Pi — P2 = fd + dp where r? is a strongly contact (n + A;)-form 
(i.e., such that pkfd = 0) and p is a strongly contact {n + k — l)-form {pk-iP = 0). Thus 
PkPi = PkP2 + PkdpkP PkP2- However, pi - p2 G 0 ])+^, so that 

X(pi)=X(p 2 ). 

Hence, every class [p] G is completely determined by a unique canonical source form. 

If we denote by XA))_|_^ the image of A])_|_^ by X, we can see that the interior Euler operator 

X: A;+, ^ XAO+fc C p ^ X(p), 

induces a representation mapping, 

Rn+k-. v:+k^K7k^ R„+fc([p])=x(p). 

Eor k = 1 there is no other source form representing a class [p] G so that every class is 

completely determined by a unique dynamical form. If [p] = £n{p) = [dp\ thenX(fip) = I{dhp) is 
the Euler-Lagrange form of the Lagrangian hp. Note that this argument proves the uniqueness 
of the Euler-Lagrange form. At the same time Proposition 3.1 shows us that for k > 2 the 
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representing canonical source form is no longer a unique source form in the class [p]. An explicit 
important example of nonuniquness for fc = 2 is realized by Helmholtz forms, and will be 
discussed below in Section 5.1. 

• If g = 0 we define Rq be the identity mapping. 

Now, by means of the representation mappings Rq we construct to the variational sequence 
0 —)■ My —)• a representation sequence 

0 ^ My ^ Ro(Vo") ^ Ri(V[) ^ R2{VI) ^ • 

shortly denoted by 0 —)■ My —)■ and called the Tokens representation of 0 —)■ My — )■ V^. 

We also denote by 


Eq: R,(V;)^R,+i(V,V), q>l, 

the morphisms in the representation sequence. Hence, the definition of the Eq’s follows from 
the commutativity of the diagrams 


£a- 


9-1 , 


yr 

<? 


yr 

'^q+l 


-9 + 1 


Rn 


Rq + 1 


(3.4) 




Rq+liV^q+,] 


^9+1 , 


and we can immediately see that the following theorem holds true: 
Theorem 3.4. The representation sequence 0 —)• My —)• is exact. 

Indeed, by definition of Eq and Sq, 


Eq ° Rq{[p\) = Rq+l ° ^qiip]) = Rq+l{[dp])^ 
and for every p £ Rq{Vq) we have p = Rq{[p]) for some [p\ G Vq. Hence 

{Eq+l O Eq){p) = [Eq+l O Eq O Rq){[p]) = (Eq+l O Rq+l O £:,)([/>]) 

= {Eq+l O Rq+l){[dp]) = {Rq+2 O Sq+li[dp]) = Rq+2i[ddp]) = 0. 


The representation sequence can be easily understood, and explicit formulas for the mor¬ 
phisms Eq and the forms Eq{p) G Rq-\-i{Vq^i) can be easily obtained, if we write the diag¬ 
ram (3.4) in the explicit form 


A^/er a ;;+ 2 / 0;;+2 —^ 


Rn 


1 


Rn + l 


Al'+i 


^n +1 


Rn + 2 


1KV2 


where Rq, 1 < q < n, has the meaning of horizontalization h on classes of order r, and Rn+k is 
the operator I acting on pi^p where p is of order r. Then since elements of the sheaves Rq{Vq) 
are functions for q = 0, horizontal forms for 1 < g < n, and canonical source forms for q > n-|-1, 
we have 


Eoif) = hdf, Eq{hp) = hdhp, 1 < ^ < n — 1, 
En{hp) = Z{dhp), En+k{R{p)) = R{d{T{p))), k>l. 
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In particnlar, for q = n we can write in coordinates hp = Lloq, where coq = dx^ A • • • A dx^, 
and then 


En{hp) = Rn+i{[dp]) =d:{dp) =I{dhp) = Y, i-'^ f^dj 

\J\=0 



A wq. 


Hence, En- Rn{V^) —)■ -Rn+i(V^_,_;^) is the Euler-Lagrange mapping, assigning to every Lag- 
rangian X = hp its Euler-Lagrange form Ex = X{dX). The next morphism, E^+i : Rn+iiVJ^j^i) —)• 
i?ri+ 2 (V^+ 2 ) Helmholtz mapping, assigning to every dynamical form e its canonical Helm¬ 

holtz form = I{de). Note that the preceding morphism to the Euler-Lagrange morphism, 
En-i : Rn-i{V^_i) — >■ RniVJ^}, assigns to every horizontal (n—l)-form (p (resp., if dimX = n = 1 , 
to a function /) a Lagrangian A = hdip (resp. A = hdf). This is so-called null-Lagrangian (also 
called variationally trivial Lagrangian), due to the fact that its Euler-Lagrange form is by 
exactness of the sequence equal to zero. 

More generally, for q > 1, elements of RqiVq) belonging to the kernel of the variational 
morphism Eq are called variationally trivial. 

We say that two elements rji, p 2 of RqiVq) are (variationally) equivalent if their difference 
is variationally trivial, i.e., ri 2 — pi G Kei Eq. Thus, two Lagrangians are equivalent iff they 
differ by a null Lagrangian, hdp, two dynamical forms are equivalent iff they differ by a locally 
variational form, and two Helmholtz-like forms are equivalent iff they differ by a Helmholtz 
form. Generally, two canonical source forms are equivalent iff they differ by X(d{Xp)). 

Example 3.5. As an illustration of the use of the Takens representation of the variational 
sequence let us compute the Euler-Lagrange operator with help of the operator X. Indeed, given 
a Lagrangian A, the Euler-Lagrange form of A arises either as the image of A by the morphism En, 
or, by application of X to the exterior derivative of A (remind the formula En{X) = X(dX) coming 
from the commutativity relation En oX = X o d). 

Consider the Lagrangian of a free quantum particle moving in one dimension. The quantum 
mechanical state of such a system is, in the coordinate representation, described by a wave 
function (j)(t, x): M x M —)• C, and the Lagrangian is 

A = ‘^ 0 , uJo = dtA dx. 


Denoting v(t,x) = Re4>it,x) and w(t,x) = lm(j)(t,x) we obtain 


A 



h 

2 


(vwt — WVt) 


wo, 


(3.5) 


which is a first-order Lagrangian for the fibred manifold vr: x —)• (hence n = 2, m = 2, 

r = 1 ) where vr is the canonical projection, with fibred coordinates (x^, x^, y^,y‘^,y\,y\,y‘l,y 2 ) = 
(t,x,v,w,vt,Vx,wt,Wx). In this notation, the generating contact I-forms on J^(M^ x M^) read 
as follows: = dv — vtdt — Vxdx, uP' = dw — wtdt — Wxdx, uj} = dvt — vudt — vtxdx, oj], = 

dvx — Vxtdt — Vxxdx, = dvjt — wudt — wtxdx, u)^ = dwx — Wxtdt — Wxxdx. Now, pidX is the 
following 1-contact 3-form on J^(]R^ x M^): 


pidX = 


h , h , h , h ^ , 

--wtdv -h -vtdw + -wdvt - -vdwt - —Vxdvx - —WxdWx ) wq. 
2 2 2 2 2m 2m 


Calculating X(dA) using formula (3.1) we obtain 


X(dX) = A 


h 


--wt -dt{ -w] -dx{ 


h^ 


Wo 
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+ UJ^ ^ 


'-Vt - dt ( -'-V ] -d^( 


UJQ 


\ 1 / \ 2 
-nwt + I ^ + I nvt + I ^ 


A Wo- 


The dynamical 3-form X{d\) is the Euler-Lagrange form of (3.5). This means that 

ei = —Vxx - dwt, S2 = —Wxx + nvt, 

2m 2m 

are the Euler-Lagrange expressions of the Lagrangian (3.5), and equations ei o = 0, 
62 o = 0 are the Euler-Lagrange equations for extremals 7 of (3.5). Denoting e = ei + 162 

we obtain the Euler-Lagrange equation 

^2 \ 

— + ihcpt j O J^7 = 0, 


which, indeed, is the Schrodinger equation. 


3.3 Lepage n-forms and the first variation formnla 

Lepage n-forms were introduced by Krupka in 1973 [46] (see also [48]) in order to establish 
foundations of the higher-order calculus of variations in jet bundles. They are fundamental for 
a geometric formulation of the intrinsic first variation formula, and of coordinate free, global 
Lagrangian and Hamiltonian mechanics and field theory. Combined with the concepts of inva¬ 
riant variational functionals they provide geometric formulations of Noether theorems, as well 
as geometric integration methods based on symmetries. 

Here we remind only some basic properties of Lepage n-forms, explored in the variational 
sequence theory. For more details and applications we refer the reader to the survey papers 
[54, 59, 60, 67, 76], and the book [55]. 

As above, vr: V —)• X is a fibred manifold, n = dimX, and m = dim V — n. 

Definition 3.6. Let r > 0. A n-form p on J'^Y is called a Lepage n-form of order r if for every 
TTr.o-vertical vector field ^ on ,YY 

HC\dp) = 0 . 

Since p is an n-form on J'^Y its horizontal component hp is a Lagrangian on 
Notice that 

• every n-form on V is a Lepage n-form. The corresponding Lagrangian is then defined 
on J^Y, and it is a polynomial of degree n in the first derivatives^, 

• every closed n-form on J^Y is a Lepage n-form. The corresponding Lagrangian is a null- 
Lagrangian (giving rise to the zero Euler-Lagrange form). 

The structure of Lepage n-forms is characterized as follows: 

Theorem 3.7 (Krupka [46]). The following conditions are equivalent: 

(1) p is a Lepage n-form of order r, 

(2) pidp is a dynamical form, i.e., pidp = I{dp), 

^This fact is of particular importance in Hamiltonian mechanics and field theory, or in dealing with variational 
forces and energy-momentum tensors [38, 54, 68]. 
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(3) ^dp = E + F, where E is a dynamical form, and F is at least 2-contact, 

(4) = 6 hp + du + p = 6hp-\-pidiy-\-ri, where Ohp is the Cartan form of the Lagrangian hp, 
v is a contact (n — l)-form, and p, resp. p is at least 2-contact. 

It should be stressed that for n > 2 and r > 3 the Cartan form is generically not global. 
However, the above theorem states that it can be “globalized” by adding pidu. Indeed, in the 
decomposition of ^T*_^_l^p in (4), O^p -\-pidv and p are global (being the at most 1-contact and 
the at least 2 -contact part of p, respectively), while, in general, the decomposition 6 hp + pidv is 
not invariant under changes of fibred coordinates. 

The name ‘Cartan form’ refers to E. Cartan, who introduced the form to the classical calculus 
of variations [10]. Concerning different aspects of generalization to many independent variables 
we refer, e.g., to [ 6 , 9, 12, 18, 29, 33, 34, 46, 47, 78, 90]. Higher-order Cartan forms have been 
introduced and studied by many authors; see, e.g., [12, 13, 21, 22, 30, 40, 48] to name just a few. 

Accounting the definition and properties of the operators X and TZ we are able to find an 
intrinsic formula for the Cartan form: 

Theorem 3.8. Let r > 0, and p G A”. Then the Cartan form of the Lagrangian X = hp takes 
the form 

6 x = X-piTZ{dX). (3.6) 

Notice that, indeed, 0hp is a Lepage n-form, equivalent with p (and hp), since by formulas (3.2) 
and (3.3) 

VidOhp = d^idOhp) + pidpiTZ{d6hp) = 1{d9hp) + pidpiTl{dhp) - pidpiTZ{dpiTZ{dhp)) 

= l{d6hp) + pidpiTZ{dhp) - pidpiTZ{dhp) + T{dpiTZ{dhp)) =X{d6 hp), 

and, writing ^p = hp -\- j3 we can see that (for a proper s) 

K,rP - dhp = l5+piR{dhp) G 0^. 

The meaning of Lepage n-forms for the calculus of variations comes from the fact that if p 
is a Lepage n-form then the dynamical form pidp = E^p is the Euler-Lagrange form of the 
Lagrangian hp. Indeed, the components of pidp are the Euler-Lagrange expressions of the 
Lagrangian A = hp. 

We stress that even though Ohp need not be global, the local forms do give rise to the global 
form pidOhp. Indeed, formula (3.6) yields 

PidOx = pidX — pidpiR,{dX) = Z{dX), 

which, indeed, is a global form. (Note that we recovered the formula for the Euler-Lagrange 
form of A in Takens representation of the variational sequence.) Moreover, since for every Lepage 
form p, rPidp = pidOhp, we get another proof of uniqueness of the Euler-Lagrange form. 

As the horizontal parts of Lepage n-forms are Lagrangians, we can view a Lepage n-form as 
an extension of a Lagrangian by a contact form, and we come to the concept of Lepage equivalent 
of a Lagrangian [46]. Given a Lagrangian A, a Lepage equivalent of A is a Lepage form p such 
that A = hp. 

Remarkably, if A is of order r, then its Lepage equivalents are generically of order 2r — 1, and 
the Euler-Lagrange form is of order 2r. 

Local existence of Lepage equivalents of A and their structure follows immediately from the 
above theorem. Moreover, it is well known that every Lagrangian admits a global Lepage 
equivalent. Eor n = I and any r, or for any n and r < 2 a global Lepage equivalent of A 
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is the Cartan form. If n = 1 (mechanics and higher-order mechanics) the Cartan form 9x is 
the unique Lepage equivalent of A. For n > 1, Lepage equivalent of A is no longer unique. 
Remarkably, apart from the (globalized) Cartan form, there are also other distinguished global 
Lepage equivalents of A, as, e.g., the celebrated Caratheodory form ([9, 78] for r = 1, and [92] 
for r = 2), which is invariant with respect to all (not only fibred) coordinate transformations, 
and the Krupka-Betounes form (for r = 1) [6, 47], which has the property that dp = 0 if and 
only if Ehp = 0 (i.e., the Lagrangian belongs to KerFi^). 


Example 3.9. Let us illustrate the use of the Lepage representation of the variational sequence 
on an example. Consider the quantum particle Lagrangian (3.5) in Example 3.5 (including 
notations). Its Cartan form reads 

/I T 1 1 2 1 1 1 2 7 

9\ = Lujq -|- t:—w f\dx + -—uj f\dx — ——uj Adt — -—uj A dt 

OVf UWt 

= (- ^ + wl)-^ {vwt - tuut)^ wo 

1 , 2 7^17 ^27 

H- VxOJ Adt -\ - WxOJ Adt -\—woj Adx - vu A dx. 

2m 2m 2 2 


As we have seen, being a Lepage equivalent of the Lagrangian (3.5), the form 9\ quickly gives 
us the Euler-Lagrange form E\, and hence the Euler-Lagrange equations: 


Ex = pid9x = —Vxx - hwt 
2m 


ui A Wo -|- 


—Wxx + nvt 
2m 


w A Wo. 


We can see that, indeed. Ex equals to I{dX) obtained in Example 3.5; this follows from the 
commutativity relation Lep 2 0 E 2 = do Lep^. 

We note that the form 9x (or any other Lepage equivalent of the Lagrangian) also gives 
local Noetherian conserved currents related with symmetries of the Lagrangian and the Euler- 
Lagrange form (see the discussion below). 


With help of a Lepage equivalent of a Lagrangian the non-intrinsic procedure of integration 
by parts in the first variation of the action function is substituted by a geometric splitting 
with help of the Cartan formula for decomposition of the Lie derivative. As a consequence, the 
integral first variation formula appears also in the following intrinsic differential form [46] (for 
more details see, e.g., [54]). Given a Lagrangian A of order r, and a vr-projectable vector field H 
on Y (a “variation”), it holds 

Ljr^X = Lj2r^hp = h(j2”-lHjdp) + /ld(j2”-lHjp) = hLj2r-l^p, (3.7) 


where denotes the Lie derivative along the r-jet prolongation of H, and p is (any) Lepage 
equivalent of A. Let us show that the first term in the sum of (3.7) is the Euler-Lagrange term 
(carrying information about extremals), and the second term is the Noether term, carrying 
information about conservations laws. 

We have A = hp, and the Euler-Lagrange form is defined by Ex = pidg. For the first 
differential form on the right-hand side of (3.7) we obtain 

h(j2”-iHjdp) = h(j2”Hj7r2%2r-i^F) = h{,f^E\{pidp + p)) = h{J^^E\iEx)), 

where p is at least 2-contact, hence is contact and /i(J^”Sjp) = 0. Denote Ex = 

E^jU^Acoo (he., E^j are the Euler-Lagrange expressions of X),uji = ^Jwq, and H = 

Then 


^E\dp) = /i(J^”Hj {EfjOj^ A wq)) = h^EfjE^uo — A Wj) = Eq-H'^wq, 

showing that, indeed, h{J‘^^~^E\dp) provides Euler-Lagrange equations. 
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Let us turn to the second differential form in the sum of (3.7). Assume that H is a Noether 
symmetry of the Lagrangian A, meaning that Ljr’^X = 0. Then in the infinitesimal first variation 
formula (3.7) the left-hand side vanishes, and along every extremal 7 of the Lagrangian A the 
n-form h{j'^^~^E\dp) vanishes, because the Euler-Lagrange equations Eo-o = 0, 1 < cr < m, 
hold true. Consequently, the n-form hd{j‘^^~^'B\p) vanishes along extremals, which is the cele¬ 
brated Noether theorem. The (n — l)-form 

$(5) = 


is then the corresponding Noether current. 

In the simplest case of first-order mechanics when p = 6\ 
well-known formulas 




/ (9L d dL \ 

\ dq^ dt ^ 


cI)(H) 






-He+Pa'S^. 


Ldt + we recover the 


If H is a Noether symmetry of A, then the function <h(H) is constant along extremals. 

We shall return to the variational splitting of the Lie derivative in a more general context in 
Section 4. 


3.4 Lepage forms of higher degrees, and the Lepage representation 
of the variational sequence 


A motivation for introducing Lepage forms of higher degrees is the following observation: If 
dimX = 1 then every Lagrangian A has a unique Lepage equivalent, the Cartan form 9\. 
The mapping Lep^: A^ 9 A —)■ 0 a £ thus relates the Euler-Lagrange mapping to the 

exterior derivative through the relationship E\ = pidOx- This suggests the idea to extend Lep^ 
to dynamical forms, i.e., to Lep 2 : A^ 9 e —)> Oe G A|, so that we would have 


A 


El 


■> 


Ex 


Lepi 


Lepa 


Ox dOx 


The concept of a Lepage equivalent of a dynamical form was introduced in [64] for locally 
variational forms. In that situation one obtains an extension of the Euler-Lagrange form to 
a closed form which (similarly as in the case of a Lepage equivalent of a Lagrangian) is unique for 
re = I and nonunique otherwise (see [37, 38, 64, 66 , 72]). Lepage equivalents of Euler-Lagrange 
forms became important particularly in study of the inverse variational problem, Hamiltonian 
theory, and geometric integration based on symmetries of the equations [34, 55, 65, 66 , 68 , 71]. 

Having the variational sequence and its Takens representation allows to define the concept 
of Lepage equivalent for any source form, and, in this way, to link all the variational operators 
to the exterior derivative. Remarkably, with Lepage forms one gets another representation 
of the variational sequence where variational morphisms simplify to exterior derivatives. The 
immediate benefit is the reformulation and solution of the problem of existence of Lagrangians 
for given differential equations “as they stand”, the problem of the structure of null Lagrangians, 
and the other inverse problems in the variational sequence: by using Lepage forms these problems 
are reduced to application of the Poincare lemma. 

The idea how to generalize the concept of Lepage form to higher degrees is suggested by 
property (2) (or (3)) in Theorem 3.7, that the lowest contact component of the exterior derivative 
of a Lepage form p is the canonical source form for dp: 
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Definition 3.10 ([45, 61]). Let A: > 0. A (n + /c)-form p on J^Y is called Lepage form if 


Pk+idp = I{dp). (3.8) 

Lepage (n + /c)-forms for k > 0 have many similar properties as have Lepage n-forms. First, 

• every g'-form on V, g > n, is a Lepage form, 

• every closed g-form on g > re, is a Lepage form. 

Structure of Lepage forms can be derived by direct calculations in coordinates from the 
definition. However, this procedure is quite lengthy and tedious. Much more advantageously, 
we again explore the operator TZ to solve the equation (3.8) with respect to p in an intrinsic 
way. Then we obtain: 

Theorem 3.11. Equation (3.8) has the solution 

K+i,rP = ^PkP + du + p = Op^p + Pk+idiP + r], (3.9) 

where 

dpkP = PkP - Pk+i'hZidpkP), (3.10) 

re is an arbitrary at least {k + l)-contact {n + k — l)-form, and p {resp. rj) is an arbitrary at least 
{k + 2 )-contact form. 

For every choice of v and p {resp. rf) the Lepage forms (3.9) belong to the same variational 
class [p] E {for a proper s), i.e., 7r*,.p - Op^p E 0^+^. 

Proof. First we show that 6p^p given by (3.10) is a Lepage form, i.e., satisfies Pk+id9p^p = 
X{d6p^p). Using the same procedure as above for /c = 0, we obtain: 

Pk+id0pj.p = ^{ddpi^p) + Pk+idpk+iT^{d0pi,p) 

= X{d0p^p) + pk+idpk+i'R,{dpkp) - Pk+idpk+iTZ{dpk+i'R,{dpkp)) =X(d0p^p), 


since 


Pk+idpk+i'JZ{dpk+i'JZ{dpkp)) = pk+idpk+i'R,{dpkp) -X{dpk+iTl{dpkp)) 

= pk+idpk+iR{dpkp), 

in view of X{dpk+iR.{dpkp)) = X{pk+i{dpk+iTZ{dpkp))) = 0 . Hence, 6p^p is a solution of (3.8). 
Since the equation concerns dp, it is clear that if p is a solution then also p + dre is a solution, 
where re is an arbitrary {n + k — l)-form. Moreover, equation (3.8) allows to determine only 
the lowest, {k + l)-contact component of dp. Thus, any solution p is determined up to dre + p, 
where p is at least {k + 2)-contact. In this way, so far we have obtained p = Op^p + dre + p, 
where p is at least {k + 2)-contact. We observe that Pkdp^p = PkP-, hence Pfcdre must be equal to 
zero, meaning that re is at least {k + l)-contact. 

Finally, we have to show that whatever the choice of re and p, p is equivalent with Op^p. This 
is, however, easily seen: if we denote by s the order of Q-py.p, and by f5 the at least {k + l)-contact 
component of p, we obtain tt* ^p — Op^p = pkP + /3 — 9p^p = fd -d-pk+iTZ{dpkp), which is a strongly 
contact (re + A:)-form, thus belonging to the kernel 0® _!_/.• ■ 

The (re + A:)-form 9pj^p is the higher-degree generalization of the Cartan form. It is completely 
determined by its lowest contact (/c-contact) component. Similarly as the Cartan form of degree re 
it is generically not global, since it is defined by means of the operator TZ. However, quite 
similarly as in the case of Lepage re-forms, it can be “globalized” by adding a proper termp^+idre. 
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And similarly as in the familiar k = 0 situation, Pk+idOp^p is global, and it holds 
Pk+idOp^p = I{dpkp). 

Also in the higher degree situation, there is a distinguished case of n = 1 (and arbitrary r) 
(mechanics). Then v a. /c-form and p is a {k + l)-form, hence they both are zero and we have: 

Corollary 3.12. Let n = dimAi = 1, r > 0. Then for every k > 0, the higher degree Cartan 
form Op^p is unique and hence global. 

The higher-degree generalization of Theorem 3.7 now reads as follows: 

Theorem 3.13. Let k >0. The following conditions are equivalent: 

(1) p is a Lepage (n -|- k)-form of order r, 

(2) ^dp = 1{dp) + F, where F is at least {k + 2)-contact, 

(3) pk+idTZ{pk+idp) = 0, 

(4) K-ei,rP = ^Pkp + dv + p = 9p^p + Pk+idu + rj, where Op^^p = pkP - Pk+iF.{dpkp), v is an at 
least [k -\- l)-contact {n + k — l)-form, and p, resp. p is at least {k + 2)-contact. 

Again, in the decomposition of in (4), Opj^p + pk+idu and p are global (being the at 

most (A;-|-l)-contact and the at least (A:-|-2)-contact part of p, respectively), while, in general, the 
decomposition Op^p -\-pk+idv is not invariant under changes of fibred coordinates. And pk+idp 
is unique (independent upon a choice of p). 

Also the concept of Lepage equivalent of a Lagrangian extends to (n -|- A:)-forms. Given 
a /c-contact form u, by a Lepage equivalent of a we mean a Lepage form p such that pkP = cr. 
Theorem 3.13 then gives the structure of Lepage equivalents, and guarantees local existence. 

There remains to answer a question about existence of global Lepage equivalents of higher 
degrees. The result is affirmative, and its proof is a generalization of a proof for n-forms [48]. 
The idea is to show that there exists a Lepage equivalent which is a global form on a closed 
submanifold of Then, since the (local) form defines a soft sheaf of forms it can be 

extended to a global form on the whole space. 

Let us consider the canonical injection ir-+i,r : J'’'''^(J’’T^) defined by ir-eixijx^^^p) 

=jrn^’’7)- 

Theorem 3.14. Let cf he a {n-\-k)-form such thatpicf = 0, I > k-\-l, e.g., (f is the pk component 
defined on J^{J''Y) of a {n -\- k)-form on T'Y. There exists a (n -|- k)-form p on J'^~^^{J'^Y) 
such that 

(1) p is {'Kr)r+i,o-horizontal and p is a {k l)-contact (n -|- k)-form on J''^^{J'^Y); 

(2) for every {'Kr)r+ifl-VGrtical vector field Z on {i.e., such that the projection 

of Z on J^Y is zero) 

C+l,rLz{Mr+pi^ + P) 

is a {k + l)-contact (n -|- k)-form on J''^^{J'^Y); 

(3) the form r((^r’)r-i-i i*/* F p) is a Lepage (n -|- k)-form. 

Proof. It is enough to prove uniqueness of p when restricted to the snbmanifold ir+i,r{J‘^'^^^Y). 
In fact, for pa {n-\-k) form on J''Y, we can take p = —X{pkp) when </> = pkp. The result follows 
by the uniqueness of X{pkp) which implies that it is a globally defined form on a closed subset 
of (defined by the injection above), therefore since the (local) form defines a soft sheaf 

of forms on J'^^~^^Y, then p can be globalized. ■ 
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The most important Lepage forms of degree > n are Lepage equivalents of canonical source 
forms. By the above, if u is a canonical source form of degree n + k, k > 1, and order r, i.e., 
X((t) is 7r2r+i,r-pi’ojectable, and pkO' = a = X(it), then all its local Lepage equivalents take the 
form 


p = 0a + du + n = 6a + Pk+idv + r], 


(3.11) 


where 


0a = a - pk+iTZida), 

V is an arbitrary at least {k + l)-contact (n + fc — l)-form, and p (resp. rf) is an arbitrary at least 
{k + 2)-contact form. Moreover, 

Pk+\d0a =T{d(j). (3.12) 

In particular, if e is a dynamical form then P 2 d 6 s is the canonical Helmholtz form = I{de). 
Comparing Lepage forms of different degrees yields the following results: 

Theorem 3.15. Let q > n. If a is a Lepage equivalent of Rq{[p\) then da is a Lepage equivalent 
of Rq+i([dp\). 

Proof. Since a is a Lepage form, da is Lepage, as trivially follows from the definition. Denote 
a = Rq{[p]). q = n, a is a Lagrangian, hence R„([p]) = hp, and Rn+i{[dp]) = R{dp) = 
Z{dhp). Thus a = 0hp + dv + p, and da is a Lepage equivalent of pida. However, by (3.12), 
pida = piddhp = I{dhp) = Rn+i{[dp]). If q = n + k, k > 1, we have Rn+ki[p]) = df{p), and 
Rn+k+i{[dp]) = Z{dp) = I{dl{p)) = X{da). The form da is a Lepage equivalent of pk+ida., 
which by (3.12) is equal to pk+id6a = L{da) = Rn+k+i{[dp\). ■ 

Theorem 3.16. A closed Lepage equivalent of a canonical source form of degree n + k, k > 1, is 
locally equal to da where a is a Lepage equivalent of a canonical source form of degree n + k — 1 
if k >2, resp. of a Lagrangian if k = 1. 

Proof. If /3 is a closed Lepage equivalent of fi = pkP = X(/3), then locally (3 = da where (with 
help of the contact homotopy operator), a = Aj3, and Pk-ict = Apk(3 = Aa. 

If k = 1 then pga = ha is a Lagrangian, and pida = pij3 = X(/3), so that a = Aj3 is a Lepage 
equivalent of ha = Aa. 

If k > 2, set e = X{Aa). e is a canonical source form of degree n + k — I representing the 
class [Aa]. Therefore (by the preceding theorem), if a' is a Lepage equivalent of X{Aa) then da' 
is a Lepage equivalent of X{de) = X{dAa) = X{dpk-ia) = X{da) = X{j3) = a. Now, both da' 
and /3 are Lepage equivalents of a, hence we have f3 = da' + dv + p, where u is at least {k + 1)- 
contact and p is at least {k + 2)-contact, and since dj3 = 0, we get dp = 0, i.e., locally p = dr 
and T = Ap. Hence, ft = d{a' + v + t), and a = M/3 = a'+ i^ + r is a Lepage form (differing from 
the Lepage {k — l)-form a' by an at least [k + l)-contact form), and it is a Lepage equivalent 
of Pk-iot = Pk-iOi' = X[Aa), as required. Moreover, we note that since Pk-iOi = Apkfl = Aa, it 
holds X{Aa) = Aa = A{Xa). ■ 

Corollary 3.17. 

(1) For every canonical source form a 
X{Aa) = Aa = A(Xa). 


(2) If a is a canonical source (n + k)-form then Aa is a canonical source {n + k — l)-form. 
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Corollary 3.18. 

(1) If (5 is a closed Lepage equivalent of a dynamical form e then A/3 is a Lepage equivalent of 
the Lagrangian A = Ae. 

(2) If (3 is a closed Lepage equivalent of a canonical source form a then A/3 is a Lepage 
equivalent of the canonical source form Aa. 

The above properties of Lepage forms suggest a new representation of the variational se¬ 
quence, based on Lepage equivalents of canonical source forms, as follows: For every k > 0 we 
have a mapping Lep„_,_;. assigning to every Lagrangian A G , if /c = 0, and every canonical 
source form a G X(A^_|_^) if A: > 1, the family {p} of its Lepage equivalents (3.11). 

With help of Lepage mappings Theorem 3.15 can be reformulated as follows: 

Theorem 3.15'. If Lepq{Rq{[p])) = {a} then Lepg_,_]^(i?q+i([d/?])) = {da}, or 

{dLepqiRqi[p]))} = {Lepg+i(R5+i([d/9]))}. 

The Lepage maps then induce the following representation mappings: 

Rq = Rq for 0 < g < n — 1 and Rq = Lep^ oRq for q > n. 

Note that for g > n the images by Rq of classes in Vq are subsheaves of A^. 

With help of the representation mappings Rq we construct to the variational sequence 0 —)• 
My —>■ a representation sequence 

0 —)• My —)• Rq{Vq) —)• Ri(V]') —)• R2iV2) ■ 

shortly denoted by 0 —>■ My —)• R^/V}), and called the Lepage representation of the variational 
sequence 0 —)■ My —)• V{. 

Theorem 3.19. The representation sequence 0 —>■ My —)• R*(V^) is exact, and every morphism 
^qC^q) ^q+iO^q+i) q > n, is the exterior derivative operator d [acting on classes of 
Lepage forms). 

Proof. Denote 

Eq: i?,(V;)^R,+i(V,V), g>n. 

With the representation mapping Rq, Theorem 3.15 states that ‘if Rq[[p]) = {a} then Rq^i[[dp]) 
= {da}\ in other words, the exterior derivative operator d extends in an obvious way to classes 
of Lepage forms; denoting the exterior derivative of classes by the same symbol d, we can see 
that 


Eq{{a}) = {da} = d{a}. 

This proves that Eg = d, q > n. 

In view of this result, and since Rq = Rq for 0 < g < n — 1, and the Takens representation 
sequence is exact, it remains to check the exactness in the arrow En-i. To this end, let us write 
the variational sequence and its Lepage representation in form of the following diagram: 

'n 1 1 £n + l ^ 

^ '^n+1 ^ ■ ■ ■ 

Rn + 1 

^ Rn+l{Vf^,) ••• 


£n-2 


E71.-2 


■> V'’ 1 

'^n—1 
En — l 


Ki-i,x 


£n-l 




yr 

»' n 


Rn 


> Rn[v:) 
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Then for every r] G Rn-i{V^_i) = X obtain 

{d o En-i){ri) = (do En-i)iRn-i{[p])) = {do En-i o Rn-i){[p]) 

= {doRnO £n-i){[p]) = d{Rn{[dp])) = d'Le^^{Rn{[dp]) 

= ljeY>n+i{Rn+i{[ddp]) = 0 . 


It is worth to write all the three sequences within one scheme as follows: 


En-l 


> a;/©; 

Rn 


En-1 


-)• 


A^+i 


Lep„ 


Rn-\-l 

_v A^+1 

Lep^+i 


-n+l 


-n + 2 


> a;;+2 /0;;+2 —^ 

Rn + 2 


En + 1 


-^^^n+2 


En+2 


Lep, 


'n + 2 


-Bn.-l 


* Lep(A;«) Lep(A;« ,,,) Lep(IA;A) ... 

As we have seen, for n = dim A = 1 (mechanics), Lepage equivalent is unique for every k >0 
(recall that this is the Cartan form of degree n + k). This means that for every k, 


Lep„+fc : Z{h*^j^y.) —^ Ki+k^ 

i.e., Lep„_,_;i.(Xp) is an element of the sheaf A*^^ (rather than a family of elements of the sheaf), 
and d is the “true” exterior derivative of differential forms. Hence “morally” (up to orders of the 
sheaves) the Lepage representation sequence becomes a subsequence of the de Rham sequence. 

Finally we want to stress that in view of the two corollaries of Theorem 3.15 the contact 
homotopy operator A restricts to canonical source forms. This means that one can obtain 
primitives (inverse images) not only of closed Lepage forms in the Lepage representation, but 
also of canonical source forms in the Takens representation. For a primitive of a variational 
dynamical form e (i.e., a Lagrangian) one has A = Ae which in coordinates is the celebrated 
Tonti Lagrangian [94, 95]. An analogous formula holds then, indeed, also for forms of higher 
degrees: if u = En+kid) then Aa = rj. 

Example 3.20. Let us write down explicit examples of Cartan forms of higher degrees in me¬ 
chanics (n = 1). Denote local fibred coordinates on Y by {t, q^) and the associated coordinates 
on J^Y by {t, q'^, q^, q^, and let 

uj^ = dq^ - (fdt, Co^ = ojI = d(f - (fdt, Cb^ = uj^ = dif - q'^dt. 


By Theorem 3.11, formula (3.10) we have: 

• For p of degree n = 1 and such that hp = X = Ldt is of order r = 1: TZ{dX) 
hence we obtain 


0 a = A + 


dL 

dq‘ 


-UJ 


dL , ,c 


which is the classical Cartan 1-form. 

• For p of degree n -|- 1 = 2 such that pip is a second-order dynamical form, i.e., pip = e = 
EijLO^ a dt, the Cartan 2-form is 0^ = e — TZ{de). To compute TZ{de) we need the following 
formula for the residual operator applied to a second-order 2-contact 3-form: denote 


E 

i,j=0 


AUj A dt, 


AV — _ AJI 


a = 
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Then 

K(a) = ^ Aii, - A£ - I {Aji, - AS ,)) a." A a." 

-1 (AS„ - ASM A M + I {Aii - ASM A ^ 

(summation over j = 0,1, 2). Taking a = de we obtain 


I 




d dE, 


2 V dcE dt dq'' 


BE 

uj"" Nijf - ^0;“^ A lE' 


dq’' 


For a locally variational form e this is the Lepage equivalent of e, introduced in [64] (see also [65]). 
• For p of degree n + 2 = 3 such that P 2 P is a Helmholtz-like form denote P 2 P = V — 
A A dt -|- Auj'^ a dt + Auj'^ A dt; the Cartan 3-form is 0^ = r/ — TZ{dp). 

The formula for the residual operator applied to a general second-order 3-contact 4-form 


2 

«= E ^j^>*"Aa;jAn;^Adt, = A^, 

i,j,k=0 

reads 

K(a) = - ({Aii*} - slTAl) w" AwJ A< + 2{Ai?,}c.,‘’ Aw" 

- {ASSP<^i Auj^ Au4, 

where {A]^!^p} = ^{A^Svp + Apl^i, + Ai^pa), and summation over j, A: = 0,1, 2 applies. Now we take 
a = dp, so that 




d^\ 

dqP dq^ j ’ 



and {A^'Jup} 


0 for i,j, k ^ 0, and obtain the Cartan 3-form of p as follows: 


e, = T( - mdS) = I? + ([{A™} - Ui” Au" A U)* 

+ 2 ([{Aii}} - {Alfp + t{A“i} - 4{A“i}) a,” A .A- A 

+ 2 ([{Aiii} - {AfJP - }) A," AW-AW'’- 2{ AjiiJlA” A A,- A w- 

+ 2 ({ A“i} - {AlSP) w-AW-AW-+ 2 ({ ASii) - {A“i}) a” A a- A w- 


with the non-zero H’s as above (cf. [69] where a corresponding formula for the variationally 
trivial case was obtained). 


3.5 The variational order 

A key issue of the Krupka variational sequence for the variational calculus, making the main 
difference between his and other approaches (C-spectral sequences, the variational bicomplex 
and infinite order variational sequences) is that the variational sequence ( 2 . 2 ) fixes (and thus 
defines) the order of variational problems in the sense of the following definition: 
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Definition 3.21. Let r/ G be a source form, or a Lepage form of degree q > 1 defined 
on (an open subset of) J^Y. We say that rj has the variational order r if it comes from the 
variational sequence of order r (i.e., it is a Takens or Lepage representation of a variational class 

[p] e 

How to understand this concept? A variational class [p] G Ag/©^ concerns local differential 
forms of degree q and order r. Takens representation (and thus also Lepage representation) is 
based on the horizontalization operator h = pq and the interior Euler operator I which provide 
a differential form of generally higher order. For example, if r = 1 we have (for q = n) the 
Lagrangian A = Rn{[p\) = hp, (for q = n + 1) the dynamical form e = Rn+i{[p]) = T(/j), (for 
q = n + 2) the Helmholtz-like form H = Rn^ 2 {[p]) = 2i(/9), etc., all of order 2. As given by 
the structure of X, generically, the components of the differential (n -|- A:)-forms {k > 1) arising 
from AO_|_^/0))_|_^ are polynomials in the derivatives starting from the order r -|- 1 up to the 
order 2r in the linear term. Notice that such a polynomial behaviour is known to be typical 
for functions known as “variational derivatives”. However the point is that the “variational 
derivatives” are not just general polynomials: the analysis of X shows that the coefficients of 
these polynomials have a determined structure of derivatives of the components of p (which are 
of order r), with prescribed symmetrization/skew-symmetrization rules in the indices, in other 
words they have to satisfy certain symmetry and integrability conditions (“order reducibility 
conditions”). For instance, a dynamical form e of order 2 (which defines a system of seeond 
order differential equations)"^ can be of variational order 1 or 2. In particular, if e is locally 
variational (comes locally from a Lagrangian, i.e., its image by the Helmholtz morphism in the 
Takens representation sequence vanishes) we obtain that if the variational order of e is 1, we 
have local Lagrangians arising as X = hp from first order p; thus A is of order 2, or of order 1 
if hp is 7r2q-projectable. If the variational order of e is nontrivially 2, we have local 3rd order 
Lagrangians arising as X = hp from second order p, or second order Lagrangians if hp is 7r3^2- 
projectable (which, indeed, cannot be reduced to the first order by extracting a variationally 
trivial Lagrangian). The same order discussion concerns global Lagrangians if = {0}. 

Example 3.22. Lagrangians, dynamical forms and Helmholtz-like forms of variational order 
one: 

• Let n = dim AC = 1. Any p G Aj/0j takes the form 
p = Adt -|- + Cadeff ~ Adt + C^dq^, 

so that 


X = hp = Ldt, where L = A + Ca<ff ■ 


We can see that a generic Lagrangian of variational order one in mechanics is of second order, 
affine in the second derivatives. Then, in particular, reducible second-order Lagrangians (and, 
of course, first-order Lagrangians) correspond to those classes [p] which are generated by 'Kip- 
horizontal forms. 

If n = dimAC > 1, we denote wq = d/dx^^\uo, etc.; we have 

n 

P^Auo + Y, A • • • A dyy A cuq...,,, 

k=l 


where the coefficients 0^,','^^'"^'° are totally antisymmetric in the indices R .. .ik, and totally 
symmetric in the pairs (c’'i,ji)... {crk,jk)- Now, X = hp = Lujq, where 


L = A + yZ+... + Ci 


-i ^jlh ^jnin 


^In the Takens representation e is an element of the sheaf of differential forms 
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(here the C’s are just corresponding multipliers of the C’s). Thus, in field theory, a generic 
Lagrangian of variational order one is of second order, polynomial of degree n in the second 
derivatives, with coefficients obeying the symmetry-antisymmetry relations described above. In 
particular, reducible (and first-order) Lagrangians correspond to generating forms p that are 
TTi^o-horizontal. The most famous Lagrangian of this kind is the scalar curvature, giving rise to 
the Einstein equations. 

• Let again n = dimX = 1. Any p G A 2/02 takes the form 

p ~ f\dt + A dq'^ + A dq^ 

~ A {A„dt + B„iydq‘') + C^udcf A dq'^ , 

where = —Cya- In Takens representation we get the dynamical form X{p) = e = A dt, 

where 

d - BC 

E^ = A^ + B^yif - ^(2^..^) = Al. + B^yif - - 2C^yqr 

Thus a generic dynamical form of variational order one in mechanics is of order three, affine in 
the third derivatives, and polynomial of degree two in the second derivatives {with a special form 
of coefficients as above). Second-order dynamical forms of variational order one correspond 
to generating forms p such that C^y = 0, meaning that the corresponding variational class 
is represented by a uj^-generated form p and the differential equations are affine in second 
derivatives: 


p ~ Ao-cv^ A dt -h BfjyUJ^ A dq‘', i.e., E^r = A^, -h Bfjyif. 

In particular, if variational {[dp] = 0, i.e., \p\ = [dX]), they come from second-order Lagrangians 
of the form 


L — A Ccrif, 


dc„ _ dCy 

dqU Qqcr ’ 


that, indeed, are reducible (equivalent with first-order Lagrangians). 

Einally, all first-order dynamical forms are of variational order one and they arise from the 
class [p] = [AaOJ^ A df\. 

• Any p G Ag/Bg takes the form 

p ~ Auj'^ A dt + A?^yijJ^ A doj’^ -\- A^,^dq^ A doj'^ 

-\- Aoj'^ A dq^ + A dcf' A dq^ -\- Caypdcf A dq^ A dq^, 

where = -Af^, = -Al^, Bf^^p = -Bl^p, B^^^p = -B^p,,, and the C’s are totally 

antisymmetric. Then the corresponding Helmholtz-like form becomes T(p) = p = Aio’^ A 

dt + Au)'^ A dt + Au'^ A dt, where, as expected, the components have a polynomial 

structure in second and higher derivatives: 




Hi 


= Al 


1 

2 


+ BlypijP + - Al, - 2{Bl,p - Bl,p)qp) - \^{Aly - SC^ypqP) 

+ + OilupXQ^Q^ + OltupXKQ^q^q'^ + Oiaup ^ ^ + P^aupXQ^ 9 ^ + -C^ypP^, 

{Al, + Al,) + {Bl,p + Bl^p) qP, Hi, = AL - 2,C^upq^. 


The requirement that the class [p] is a;°'-generated gives A^, = 0 and C^yp = 0, and simplifies 
the Helmholtz-like forms to H^, = 0, and 


H^ =a^ 


a, 


aup 


Q + Oi] 


GVp. 




< 7 Vp 


qP. 
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If, moreover, is a Helmholtz form (i.e., \dp\ = 0), related with a w^^-generated class then it 
corresponds to a second-order dynamical form e = EaUj^ A dt, where + Bauq'^, and 

Ecru — Byfj. 

The above discussion illustrates that the range of applications of the variational sequence is 
wider than that of corresponding infinite order constructions or even of some other known finite 
order constructions. Remarkably, it gives solution of the 

Inverse order problem (order reduction problem): If a differential q-form rj of order s is 
{locally/globally) trivial, or {locally/globally) variational {its image under the corresponding mor¬ 
phism Eq vanishes) what is the ‘minimal order’ of the {local/global) preimage {q — l)-form p 
such that T] = Eq-i{p) ? 

Note that this question is answered by the variational sequence, but is not answered, e.g., by 
the variational bicomplex, or even by Takens representation(!), since p is an element of a sheaf 
of differential forms of order s which contains forms of different variational orders. 

Example 3.23. Consider a dynamical form £ (nontrivially) of order three over a fibred manifold 
with n = dimX = 1. e belongs to and represents a system of 3rd order ordinary differential 
equations. Takens sequence gives us that if locally variational, e comes from a local Lagrangian 
of order 3. However, from the variational sequence we can get information about the variational 
order of e, and hence a finer result about the order of the local Lagrangians. Namely, the 
following cases may arise: 

• The variational order of e is (nontrivially) three. This means that in the variational 
sequence e is represented by a class [p] G A^/©!. If e is locally variational, any minimal 
order Lagrangian is of order three (not further reducible). 

• The variational order of e is two. This means that in the variational sequence e is rep¬ 
resented by a class [p] G A^/©!- If locally variational, it has a (local) second-order La¬ 
grangian. 

• The variational order of e is one. This means that in the variational sequence e is rep¬ 
resented by a class [p] G K\/Q\. As follows from the previous example, if e is locally 
variational then the Tonti Lagrangian (which is of order three) is locally reducible to 
a second-order Lagrangian affine in the second derivatives. 


4 Lie derivative in the variational sequence 

One of the most important features of the geometric formulation of the calculus of variations in 
jet bundles is the fact that variations can be described by Lie derivatives of differential forms 
with respect to prolongations of projeatable vector fields. A famous example is the first variation 
formula which takes the integral form [46] (cf. [32]) 

[ r-f*Ljr^X= f 

where p is a Lepage equivalent of A, H is a projectable vector field on Y (“variation vector 
field”), and H C lA is a compact connected n-dimensional submanifold with boundary. It is, 
however, even more interesting (and fundamental) that with use of the differential calculus in 
jet bundles, this formula can be (as mentioned in Section 3.3), invariantly expressed in the 
“infinitesimal form” as a formula for differential forms on the manifold J^^Y [46] 

Ljr^X = h{j‘^’'-^E\dp) + hd{j‘^’'-^E\p), 


(4.1) 
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and (with the account of A = hp) it yields a remarkable property of the Lie derivative in the 
calculus of variations: 


Lj2r^hp = hLj2r-l’^p. (4-2) 

Moreover, as discovered in [62], the Lie derivative commutes with the Euler-Lagrange opera¬ 
tor, as the Lie derivative of the Euler-Lagrange form of A is the Euler-Lagrange form of the 
transformed Lagrangian: 

Lj2v^Ex = (4.3) 

Such properties of the Lie derivative provide an elegant intrinsic formulation (and an easy proof) 
of Noether theorem, Noether-Bessel-Hagen theorem [5, 77], conservation laws, and relationship 
between symmetries of a Lagrangian and its Euler-Lagrange form (for any order and any number 
of independent variables) [48, 55]. 

Having in mind the variational sequence and its representations, there arises a question about 
generalization of these properties to variational forms of any degree. 

Eirst we recall results obtained in [56, 57, 58]. Let Cr denote the contact ideal of order r. 

Definition 4.1. A vector field Z on J^Y is called a contact symmetry if it is a symmetry of 

the contact ideal Cr (meaning that for every contact form ui the Lie derivative Lzco is a contact 

form. 

Contact symmetries were characterized in [56]. It is worth notice that for any projectable 
vector field H on the r-jet prolongation Z = is a contact symmetry, and conversely, if 
a contact symmetry Z on ,EY is projectable onto X then Z = for some H on Y. 

Let k > 0. As an immediate consequence of the definition it follows that the Lie derivative 
of a /c-contact form by a contact symmetry is an at least k-contact form. This means, however, 
that the Lie derivatives by contact symmetries preserve the sheaves 0^. 

Proposition 4.2. Let Z be a contact symmetry defined on ,EY. Then for all q > 1, Lz&g C 0^. 

Equivalently we can claim that for any two forms pi, p 2 belonging to the same class [p] in 
the variational sequence, the Lie derivatives Lzpi, LzP 2 by any contact symmetry Z are also 
equivalent. Thus, following [57] we can define the Lie derivative Lz of a class 

^z[p] = [Lzp]- (4.4) 

Now one can easily prove that the property (4.3) extends for any contact symmetry, and 
through the whole variational sequence, for any morphism Sk- A^/0^ —)• A^_|_j^/0^_|_^ [57] (see 
also [58]): 

Theorem 4.3. Let a vector field Z on J^Y be a contact symmetry. Then the Lie derivative Lz 
commutes with all morphisms in the variational sequence. More explicitly, for all k > 1 

CzSkip] = £k{Cz[p]) = £k[Z\dp\. (4.5) 

Proof. Since the Lie derivative commutes with the exterior derivative, we have for any A:-form p 
on rY 

Cz[dp] = [Lzdp] = [dLzp] = [dZ\dp]. 


Writing this formula in terms of the morphism 8^ we get (4.5). 
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From now on, we shall restrict to tTs- projeatable contact symmetries (recall that then we have 
Z = for a projectable vector field H on Y). From the definition it is clear that the Lie 
derivative then preserves even individual contact components. Namely, if Z is a projectable 
contact symmetry on then Z = and 

Ljr+i^PkP = PkLjr’EP, A: > 0, 

generalizing (4.2) to variational forms of any degree. Note that we may write 

Ljr+^=.K+l,rP = Ljr+lS '^PkP = '^PkiLjrSp). 

\k=0 J k=0 

Let us turn back to formula (4.4). It means that the Lie derivative of classes of forms, i.e., 
variational Lie derivative, can be correctly defined as the equivalence class of the standard Lie 
derivative of forms, and thus represented by forms. Such a point of view opens a possibility 
to express the infinitesimal first variation formula (4.1) in terms of the variational morphisms, 
and, in this way, to obtain a generalization of the first variation formula for any degree of forms 
(with higher degree analogs of Noether theorems, indeed). 

Some results in this direction have been already achieved, namely explicit formulae for the 
quotient Lie derivative operators were provided, as well as corresponding versions of Noether 
theorems interpreted in terms of conserved currents for Lagrangians and Euler-Lagrange mor¬ 
phisms. However, only classes of forms up to degree n + 2, were considered (the latter assumed 
to be exact) [26]. The representation made use of intrinsic decomposition formulae for vertical 
morphisms due to Kolaf [42], which express geometrically the integration by part procedure, 
however, in a different way compared to that based on the splitting of the Cartan form [32, 46]. 
The decomposition formulae introduce local objects such as momenta which could be globali¬ 
zed by means of connections, and besides the usual momentum associated with a Lagrangian, 
a ‘generalized’ momentum is associated with an Euler-Lagrange type morphism. Its interpreta¬ 
tion in the calculus of variations has not yet been exhaustively exploited; as conjectured in [80], 
generalized momenta could play a role within the multisymplectic framework for field theories. 

In the sequel of this section we propose an approach to the problem of representation of the Lie 
derivative of classes of forms which encompasses previous results and generalizes to any degree 
of forms. Actually, representation of (non exact) classes of forms of any degree appears relevant, 
especially in degree n -|- 3, for the study of nonvariational differential equations arising from 
Helmholtz-like forms (related with the problem of the existence of closed n-\- 2 forms) [69, 75]. 
We shall exploit the relation between the interior Euler operator and the Cartan formula for 
the Lie derivative of differential forms. The representation of the variational Lie derivative will 
provide, in a quite simple and immediate way, ‘any degree’ generalizations of (4.1), together 
with the generalized Noether theorems, as ‘quotient Cartan formulae'. 

These results, which have an intrinsic importance from a theoretical point of view, as shown 
in Section 5, have also various concrete applications; besides the ones presented here, it is worth 
mentioning that applications of iterated variational Lie derivatives in variational problems on 
gauge-natural bundles produced physically relevant results concerning existence and globality 
of conserved quantities and on existence of Higgs fields on spinor bundles (see, e.g., [23, 24, 81, 
82, 83, 86]). 

Let H be a projectable vector field on Y. We define the interior product of a prolongation 
of H with the equivalence class of p by the formula 

.7’=JH= (4,6) 

(the equivalence class of the interior product of the vector field with the representation of the 
equivalence class of p). This definition is well posed since the representation Rq is unique. 
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Given a ( 7 -forni p defined (locally) on J’”E, we have a commutative diagram defining an 
operator Rq 

Rq{CMp\) = Rq{[Ljr^p]) = Ljs^Rq{[p]). 

This operator is uniquely defined and is equal, respectively, to the following expressions: 


Lj^^hp, 

0 < q <n, 

s = r + 1, 

Ljo-diip), 

n 1 < q < P, 

s = 2r + 1 

Lj^’EP, 

q> P +1, 

s = r, 


where P denotes the corank of the Cartan distribution on J'^Y . 

This definition enables us to deal with ordinary Lie derivatives of forms on A^, thus we can 
apply the standard Cartan formula. 

We recall a result in [43, Theorem III.11], see also [45]. 

Lemma 4.4. Let 'E be a n-vertical vector field on Y and p a differential q-form on J^Y . Then 
the following holds true for i = 1,... ,q 

J'^^‘^E\pidpip = 

and 

Ljr+2^{7rr+2,r+i)*PiP = J'"^'^E\pi+idpip + pid{J''''~^E\pip). 

As above, we denote wq = dx^ A dx"^ A • • • A dx"^. Next, let us denote 

-Ai -Ai 

It is well known that in order to obtain a representation of classes of degree re + 1 in the 
variational sequence by dynamical forms the following integration formula is used [49] 

Wjj A Wo = —dojj A uji, 

and the corresponding representation is obtained by taking the pi component obtained by ite¬ 
rated integrations by parts. In order to integrate by parts (p + A:)-forms with p < re, we need to 
generalize as 

^ A UJij. 

This enables us to generalize results given in [45]. 

r 

Lemma 4.5. Let p G 1 < p < n. Let pkP = ^ ha! Va “ 1)-contact 

|j|=o 

{p + k — 1)-forms. Then we have the decomposition 
PkP = 3(p) +PkdpkTK{p), 

where iH(p) is a local k-contact {p + k — l)-form such that 

r r 

3f(p) +pkdpk^{p) = w“ A ^ -h Yj A Ca), 

|J|=0 |7| = 1 



djha- 


with 
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Notice that djri;^ are also {k — l)-contact p-horizontal (i.e., containing the wedge of p of the 
forms dx*) {p + k — l)-forms. It is therefore well defined a local splitting of pkp, p G ^p+k^ 
1 < p < n and k > 1. 

Of course, J = I and 91 = 7^ in the case p = n. 

Example 4.6. Let p = n — 1. We can write a;“ A Cq = A w/, where some local 

/c-contact fc-forms on J^^Y. Therefore 


^ d/(u;“ A Ca) = ^ dix^^ A w; = d* ^ A w; = dj ^ d/x^^* A dx^ A LOij 

|7| = 1 |7|=1 |/|=0 |/|=0 




r—1 

j;(-l)W[^^l Ao;;,- 

A dx* = dn 

r—1 

J/|=o 


J/|=0 


r—1 

We denote by iH(p) = Y1 (—A ujij so that 
|/|=o 


^ di{uj'^ A C^) = pkdpkd\{p). 

|/|=i 

In the following we shall use such a splitting in order to split the vertical differential of 
a p-density and define a corresponding ‘momentum form’. 

Theorem 4.7. Let 0 < q < n, 'E a ir-projectahle vector field on Y and p a differential q-form 
on J' Y. We have 

1. For 0 < q < n — 1 

Rq{^J-~[p]) = Ljr+l^hp 

= W'^'^Enldnhp + ff~^^Ev\'^{dp) + dn \Pdvhp + J^~^^^H\hp), 

where a ‘generalized momentum’ is defined by 

Pdvhp = -pffR{dp). (4.7) 


2. For q = n 

Rn{Rj^E[p]) = Ljr+l'^hp = ff^‘^Ev\En{hp) + dH{j'^^^'F.v\Pdvhp + 

where 

Pdvhp = -piTi{dp). (4.8) 

Notice, that the latter formula, written as 
Ljr+i^A = J'’^‘^Ev\E\ + d_H-(/>A,s, 

where 

= J''"^^'^v\Pdvhp + ff^^EH\hp, 

is obviously the infinitesimal first variation formula (4.1). We observe that Pdvhp are momentum 
forms, and is the Noether current related with E and A. 







30 


M. Palese, O. Rossi, E. Winterroth and J. Musilova 


Proof. For all cases by the Cartan formula 
Ljr+i^hp = “J dhp + hp) 


we have 

Ljr+i^hp = {dyhp + dnhp) + [dy + dH){{J^^^'^H + J'^^^'^v')\hp') 

= d/f (hp) + J^^‘^'EH\dHhp + J'^~^'^'E,y\dvhp. 

This expression can be further characterized in more detail. 

Let 0 < g < n — 1. Since dyhp = pidp = '3{dp) + pidpiD\{dp), and since by Lemma 4.4, 

J''~^‘^Ey\pidpi^{dp) = -pod{J'""^^Ey\pi^R{dp)) = -dH{J''~^^'^y\pi'^{dp)), 

we can define a momentum form associated with a density of degree < n by (4.7) and obtain 

Ljr+i^hp = dH{r+^EH\hp) + r+^EH\dHhp + r+^Ey\{J{dp) + pidpiD\idp)) 

= r+^Ey\Jidp) + r+^EH\dHhp + dH{r+^Ey\pd^hp + r+^EH\hp). 

Let q = n. Then dnhp = 0, hence = 0. Next, again since dyhp = pidp = 

E{dp) + pidpiTZ{dp) and since by Lemma 4.4, J'^~^‘^Ey\pidpiTl{dp) = —d//(J^+^SyJpi77(dp)), 
we define a local momentum form by (4.8), and obtain 

Lj.+ishp = r+^Ey\ (X(dp) - dH{r^^Ey\pin{dp)) + dnir^^Enlhp) 

= r+^Ey\iidp) + dH{r^^Ey\pd^hp + r^^^H\hp) 

= r+^Ey\I{dp)+dHf)X,E- 

However, (fx,E is a horizontal form, and in Takens representation, X[dp) = En{hp) and dnffx^s = 
hd4>x,E = En-i{4>x,E), proving the theorem. ■ 

Theorem 4.8. Let q = n + k, k > 1. Let E be a n-vertical vector field on Y. We have 

Rn+k{k:jrE[p]) = r+^E\pk+id{i{p)) + pkd{rE\i{p)). 

Proof. For simplicity of notation, in what follows we omit the projections (so we write pkP 
instead of tt* j,_^_^pkp, etc.). 

Since Rn+k{[p]) = R{p) = PkP - PkdpkE{p), we have 

.R„+fc(^j’-s[p]) = Ljs’^Z{p) = Ljs^pkp - Ljos.Pkdpk'Eip)- 

By applying Lemma 4.4 to dpkE{p) we have 

Ljs+i^PkdpkTZip) = J''~''^E\pk+idpkdpkTl{p) + pkd{j''E\pkdpkTl{p)) ■ 

Simple manipulations show that 

Pkd{J'’E\pkdpkLi{p)) = -J^+^Sjpfc+idpfcdpfc77(p), 

so that 

Ljs+i^PkdpkUip) = 0. 

On the other hand 

Lj^^EPkP = J''^^^\Pk+idpkP + Pkd{J''E\pkp), 
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but, since pkP = '^{p) + PkdpkT^{p), by substituting we get 

J^'^^^lPk+idpkP = J^^^E\pk+id{I{p) + pkdpkTZ{p)) 

= J^^^E\pk+id{I{p)) + J^^^E\pk+id{pkdpkTZ{p)) 

= J^^^E\pk+id{I{p)) -pkd{J''E\pkdpk'R{p)). 

Again by the same substitution in pkd{J^E\pkp), we get 

Pkd{j''E\I{p)) +pkd{J''E\pkdpk'JZ{p)). 

Summing up we finally obtain 

Rn+k{CjrE[p]) = Ljs^PkP = r+^E\pk+id{i{p))+pkd{rE\i{p)). m 

By the above theorem we get a generalization of the famous Noether equation and Noether- 
Bessel-Hagen equation (i.e., equation for symmetries of a Lagrangian and of the Euler-Lagrange 
form) to a source form of any degree. Indeed, a projectable vector field S on V is a symmetry 
of Rip), i.e., Ljs'^I^p) = 0 if and only if L js+i=^X(/9) = 0, that is if and only if 

r+^Ev\pk+id{I{p))+Pkd{rEv\I{p)) = 0 . 

Theorems 4.7 and 4.8 are the desired generalization of the infinitesimal first variation for¬ 
mula (4.1) to forms of any degree. The ‘boundary term’ on the right-hand side in Theorem 4.7 
contains ‘generalized momenta’, and provides a lower-degree analog of Noether theorem in the 
same way as in the classical case q = n (k = 0). 

Lemma 4.9. Let q = n + k, k > 1, and let E be a tt- projectable vector field on Y . We have for 
every class [p\ of degree n + k 

= 0 - 

Proof. In fact, J^^^EH\dHX{p) = 0 because dHX{p) is the horizontal differential of an n- 
horizontal form. On the other hand since dvX{p) = pk+id{I{p)) and pk+idopk+id = 0, we have 
J^^^^EnldvXip) = J^~^^EH\pk-{-id(I{p)) (a (/c-|-l)-contact and (n — l)-horizontal (n-|-A:)-form); 
therefore Rn+kl^jr+iEniP]) = J''~^^X.H\pk+id{I{p))) pfcd(J®+^S hJX(p)) e 0^+^. ■ 

Since pk+id{X{p)) = X{d{X{p))) = Rqj^i{£q{[p\)), taking into account Theorem 3.15, we 
finally obtain a ‘quotient’ Cartan formula for the variational Lie derivative of n-|-A; classes, with 
k>l. 

Theorem 4.10. Let k > 1. Let E be a tt -projectable vector field on Y. For every class 
[p] G A))_|_^/0))_|_^ we have {up to pull-backs) 

Rj^eIp] = J'^'^^'Xv\£n-i-k{[p]) + £n-i-k-l{j'^'^^'Xv\ [p]) ■ 

Proof. Since H is assumed to be a vr-projectable vector field on Y, up to projections, 

Rn-{-k{Rj^E[p]) Rn+k{Rj‘r+i-'EYiP\') A Rn+kiRj^+^Ejjip]) 


holds true and the result follows from the above lemma and theorem (the latter applied for the 
vertical part of J'^E). ■ 
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5 Some applications 

5.1 Helmholtz forms in mechanics 

As mentioned above, for k > 2, a class [p] G in the variational sequence can be represented 
by different source forms. We shall illustrate this feature on the image £i„+i([e]) = [de\ G Vn+ 2 i 
where e is a second-order dynamical form, in the case when n = 1 (mechanics). 

Let n = dim A = 1; as above, denote local fibred coordinates on Y by and the 

associated coordinates on J^Y by {t, , Q'^), and let 

uj^ = dq^ - q^dt, u!^ = dq^ - (f dt, = dq^ - 'q'^dt. 

Consider a dynamical form e G A^, e = A dt. Then [de\ is represented by the following 

canonical source form (called canonical Helmholtz form), obtained by a straightforward calcula¬ 
tion: 


H,=Z{de) = - 


1 ( (dE„ dE„ 


dq'' 


+ 


dE„ ^ dE,, 


dq^ dq' 


- 2 


dq^ 
d dE, 
dt dq^ 


1 d 

2 dt 


dE„ dE, 




dcH 

BE, 


dq 

BE, 


1 d^ 

^ 2 dt 2 


BE„ BE, 


Bq^ Bq° 


Bq^ 
Aui'^ A dt. 


Bq^ 


UJ 


By construction G Ag, however, it is projectable to J'^Y. 

From the exactness of the representation sequence we get that the condition E 2 {e) = = 0 

is necessary and sufficient for existence of a local Lagrangian A such that e = Ei(A) = E\. 
Expressed in terms of the components of the condition = 0 gives the celebrated Helmholtz 
conditions of local variationality for dynamical forms, in the context of the variational sequence 
called also source constraints. 

The canonical Helmholtz form is not the unique source form providing Helmholtz conditions. 
Another distinguished Helmholtz form is [53] 


He 


2 \ \Bq^ Bq^ 2dt\ Bcf Bq^ )) ^ 

(B^ d^_ ± (BE, 

V Bq^ ^ Bcf dt V Bq’' ^ Bif ) ) ^ V d(t 



Auj^ A dt. 


He is of order 3, and is, indeed, equivalent with H^. Note that it provides Helmholtz conditions in 
an equivalent and more simple form. The relationship between the canonical and the “reduced” 
Helmholtz form is as follows [70]: 


He = H^+ P2dr], 


where 


1 d fBE„ 
1~dt 


BE, 

Bq^ 


uj’' A uj° 


It is worth note that in the most frequent case when the components of e are affine in the 
second derivatives with symmetric coefficients, i.e., 


Efj — Afj Bui/pf', B — B^, 


the canonical Helmholtz form H^ is projectable to J^Y and coincides with the “reduced' 
form H^. 
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Asking about Helmholtz forms of variational order 1, that is, coming from the variational 
sequence of order 1, we get the following result: The Helmholtz form has variational order 1 
if and only if it is projectable onto J^Y, and this is if and only if the following conditions are 
satisfied: 


dB^y _ dBgp 

dqP ~ dcH ' 


(5.1) 


Note that: 

• Conditions (5.1) are a part of the Helmholtz conditions, i.e., they are necessary for e be 
locally variational. 

• Conditions (5.1) are integrability conditions, guaranteeing existence of a function L such 
that 


QqaQqU 

(the minus sign is chosen just for convenience). L is not generally a Lagrangian for e 
(indeed, the form £ need not be variational, need not be equal to 0). However, with 
help of the Euler-Lagrange form E\ oi X = Ldt we get e = E\ — F, so that in coordinates 
equations represented by e take the form 


dL d dL 
dy^ dt dij^ 


We can conclude that the requirement that has variational order one means that the 
corresponding equations are so-called Lagrange equations of the second kind. 

• It can be easily shown that the above equations are variational, i.e., = 0 if and only 

if the force F has the form of a Lorentz-like force (as, e.g., Lorentz force, Coriolis force, 
etc.). 


5.2 Inverse problems in the variational sequence 

Perhaps the best known problems which are solved with help of variational sequences and bi¬ 
complexes are the local and global inverse problem of the calculus of variations, and the problem 
of local and global triviality of Lagrangians. Let us show how these problems appear and are 
solved with help of the variational sequence and the representation sequences. 

• Variationally trivial Lagrangians, also called null-Lagrangians are Lagrangians which give 
rise to identically zero Euler-Lagrange expressions (meaning that the Euler-Lagrange form E\ 
identically vanishes). This problem concerns the classes of degree n — 1, n and n -|- 1 and the 
corresponding morphisms £n-i '■ ^ and —)• (the Euler-Lagrange mapping). 

Consider the variational sequence 0 —)• My —>• V^, and assume that a class [p\ G satisfies 
£n{[p\) = 0. Then by exactness of the variational sequence there exists a class [rj] G such 
that £n-i{[ri\) = [p] = {drj\. The condition £n{[p\) = 0 is the local variational triviality condition 
and determines the structure of Lagrangians whose Euler-Lagrange forms vanish identically: 
every variationally trivial Lagrangian locally is a closed n-form drj modulo contact forms, in 
other words, A = hdp = dHp. In coordinates, if r/ = Puy -!-••• then L = di/L If, in addition, 
H^^Y = {0} then rj may be chosen globally defined on HY. 

In Takens representation the same result follows by the following arguments: If [p] G 
satisfies £n{[p]) =0 then A = Rn{{p\) = hp satisfies the constraint on local variational triviality 
En{hp) = 0, which means that locally hp = Rn-i[{drj\) = hdp. A corresponding local Lagrangian 







34 


M. Palese, O. Rossi, E. Winterroth and J. Musilova 


is then A = Rn{[r]]) = hr]. Condition = {0} then guarantees existence of a global rj, hence 

of a global Lagrangian hr]. 

In Lepage representation, the local triviality condition £n{[p]) = 0 means that (any) Lepage 
n-form p G [p\ (Lepage equivalent of A = hp) is closed. Hence locally p = dp, i.e., A = hdp = 
{dnp =)hdhp = dif^ujQ, and p, resp. the relevant part hp, is constructed from p with help of the 
Poincare contact homotopy operator as p = Ap, resp. hp = pUi = AX. 

• By the inverse problem of the calculus of variations, here one understands the question 
about local and global variationality of a dynamical form e. This means the problem to deter¬ 
mine constraints on variationality (Helmholtz conditions), and for a locally variational form, to 
construct a local Lagrangian. This problem concerns the classes of degree n, n -|- 1 and n -|- 2 in 
the variational sequence, and the corresponding morphisms —)• (the Euler-Lagrange 

morphism) and £n+i ■ ^n+i '^n +2 (t^® Helmholtz morphism). 

Consider the variational sequence 0 —>■ —>■ and assume that a class [p] G satisfies 

£n+i{[p]) = 0. Then by exactness of the variational sequence there exists a class [p] G such 
that £ni['n]) = [p] = [dp]. The condition £n+i{{p\) = 0 is the local variationality condition 
and determines the structure of locally variational forms (giving rise to variational equations): 
Every locally variational form e of variational order r is a closed (n + l)-form dp modulo strongly 
contact forms and exterior derivatives of contact forms of order r. Remarkably, this factorization 
procedure determines the structure of locally variational forms of the variational order r, i.e., 
produces constraints for a locally variational form of order s to come from a local Lagrangian of 
order r {order reduction constraints). If, in addition, H^^^Y = {0} then e is globally variational. 

In Takens representation, if [p] G satisfies £n+i{[p]) = 0 then the dynamical form 

£ = Rn+i([p]) = Xp satisfies E„+i(e) = 0 {Helmholtz conditions). Hence we can say: 

A dynamical form e is locally variational if and only if its Helmholtz form He = En+i{£) 
vanishes. 

Then locally e = Rn[{dp]) = Idp, providing a Lagrangian for e. 

In Lepage representation, the local variationality condition £n+i{[p]) = 0 means that (any) 
Lepage form p G [p] (Lepage equivalent of e = Xp = pidp) is closed. Hence locally p = dp, i.e., 
£ = PiP = Pidp, and p, respectively, a Lagrangian A for e is constructed from p with help of the 
Poincare contact homotopy operator as p = Ap, respectively, A = As {Tonti Lagrangian). 

In [69, 75] the same questions concerning (n -|- 2)-forms have been studied for the case of 
mechanics (dim JC = 1). Among others, explicit formulae for constraints on local ‘variational 
triviality’ of Helmholtz-like forms, £'„_|_ 2 (p) = 0, generalizing Helmholtz conditions to forms of 
degree n -|- 2, as well as Tonti-like formulas for the corresponding (nonvariational) dynamical 
forms, have been found. 

In the same way, the inverse problems, the equivalence problems, and the corresponding 
order reduction problems in the local and global form extend to each column of the variational 
sequence, and are solved by the exactness of the sequence. A practical question, however, is to 
find explicit formulae for source constraints, i.e., constraints on a source form (of any degree) 
to be variationally trivial (i.e., to belong to the kernel of the corresponding morphism in the 
variational sequence). The problem in its full generality is solved easily with help of the Lepage 
representation of the variational sequence, where it is transformed to the Poincare lemma: 

Theorem 5.1 (general inverse variational problem). Let e be a source form of degree n + k, 
k > 1. The following conditions are equivalent: 

(1) e is locally variationally trivial, i.e., belongs to KerEn+k, 

(2) £ has a closed Lepage equivalent, 

(3) every Lepage equivalent a of £ satisfies pk+ida = 0. 
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If £ is locally variationally trivial then a corresponding local primitive source (n+A: —l)-/orm rj 
{satisfying En+k-iiv) = e) is r] = As. 

If = {0} then e is globally variationally trivial, with a global primitive of order 

< 2r + 1, where r is the variational order of e. 

Proof. Let e G YeiE^+k- This means that e = T{p) where [p] G YexE^i+k- Hence [p\ = [dv\ = 
Sn+k-i{M), and we have {Lep„+fc(e)} = Rn+k{p) = Rn+k{dn) = {Lep„+fc(R„+fc(di^))} = {0}, 
proving (2). 

Next, assume d{Lepn+ki^)} = {0}- This means that every Lepage equivalent a of e satisfies 
da = dp where p is at least {k + 2)-contact. Hence pk+ida = 0, as desired. 

Finally, assume (3). Since a = 9 + du + p where p is at least {k + 2)-contact, we can see that 
Pk+id9 = 0. Hence {a} = R([/5]) where {da} = R{[dp]) = {0}. This means that [dp] = [0], i.e., 
[p] G KerT„+fc, and e = I{p) G KerEn+k- 

The remaining assertions follow immediately from the properties of A, X, and the variational 
sequence. ■ 

Condition (3) in Theorem 5.1 can be also expressed by means of the Cartan (n + A:)-form 9^ 
of the source (n + /i;)-form £ as 

Pk+id9e = 0 . 

It represents necessary and sufficient souree constraints (constraints on a source form to be 
locally variationally trivial) for source forms of any degree q > n, and of any order. For q = n + 1 
this is exactly an intrinsic form of the celebrated Helmholtz conditions mentioned above. 

5.3 Bosonic string 

A bosonic string is described as a minimal immersion of a 2-dimensional surface in a 4-dimen- 
sional pseudo-Riemannian manifold {M,g). Hence the ‘motion equation’ is the minimal surface 
equation arising from the area functional (known as Nambu-Goto action [35]). 

The corresponding fibred manifold is vr: V X where Ai = and V = X M, so that 
dim V = 6. Base coordinates and A represent the time evolution parameter and the position 
on the string, respectively, fibred coordinates (x^), p = 0,1,2,3, are local coordinates on M. 
We denote by (r*, x^, xf), where i = 0,1 and ^ = 0,1, 2, 3, associated local coordinates on AY, 
and put uj^ = dx^ — x^dH. If i: —)• S C M is an immersion of to M, and g = g^iudx^ ®dx'^ 
is a (symmetric regular pseudo-Riemannian) metric on M, we have the induced metric on S 

L*g = hijdA (g) dA = {gp,,^ o L)^:^^:^dA ® dA 


dH dA 


and the area element 


/- / / f)T^\ 

dT. = y - det{hij)dA A dA = W-det f (5^1. o i) —— j A (ir\ 

hence the Nambu-Goto action reads (in units with light speed c = 1) 


5 = -T y dS, 


where T is the (constant) string tension. In the fibred setting, we write the action as 5* = 


J^ 7 *A, where H is a domain in 


7 = 


(now considered as a horizontal form on AY) is 
A = -T det {gpiiyxfxj)dA A dA 


2 ,l) is a local section of vr, and the Lagrangian 
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(note that the action is reparametrization invariant). Let us denote 

= -TV^ = -T^J-det {g^^x'ix'<) = - gaixgpv)xQX^x{x\. 

The Cartan form of the Lagrangian A is the 2-form 

dL dL 

6\ = Ldr^ A dr^ + A dxi — A dro 

(JXq 


= TV^dr^ A dr^ + - gap9iJiu)dx^ A {x^x^x^dT^ - x^x^x^dr^). 




We introduce momenta 


dL 


^ Bxq 
dL 


pI = 


T 


Pu = 


^ dx'^ .JLId 




(gagg^i^ - ga,,g9u)XQX^x‘{^, 
{gagg^u ga/iggi^)®l2^0^0 5 


and notice that p^^Xq = p^^x^ = —T\/—D^ and = p^Xq = 0. 

As we have seen, a projectable vector field H on is a symmetry of A if H satisfies the 
Noether equation Lji^A = 0. The corresponding Noether current is then the horizontal 1-form 
'L = h{J^'B\6\), and, by Noether theorem, it is closed along critical sections (on every extremal 
surface S); d{d/ o J^y) = 0. 

To make explicit computations, denote 


d d 

= £ —r -I- -. 

5x“' 


7 ^“ — -I- ^ -I- ( d 


d^^ \ d 


^ dr'- ) dx ^' 


Contracting 9\ and taking the horizontal part we obtain 
T = S^(p°dr'-pidr°). 

This means that to every symmetry S of A we have on shell a conservation law 




dpi 

arO ^ dr^ 


= 0 . 


Symmetries of our Lagrangian can be explicitly determined from the above Noether equation. 
Since D does not depend on r® and r^, with help of momenta the Noether equation yields 


dxf^ 



We can see that the equation does not contain the projection ^ of H (which is a vector field on 
the base X = M^). This means that every vector field on is a symmetry of A, in other words, 
we recover the fact that the action is reparametrization invariant. Now se can take into account 
only vertical symmetries (which are vector fields on the space-time {M,g)). If, moreover, {M,g) 
is the Minkowski space-time, i.e., with the metric diagg = (1, —1, —1, —1), then D does not 
depend on the x^’s, and we immediately obtain that the Poincare group is a symmetry group 
of A. The symmetries and the corresponding Noether currents are (s = 1,2, 3, p = 0,1, 2,3) 


^ =-P^dT^ + pldT\ 

^ _pix°)dr° + (pgx" +p°x°)dr\ 

T = (y—p\x^ +p\x^)dT'^ d- (p?x^ — P2X^)dr^. 
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5.4 Variational Lie derivative and cohomology 

The variational Lie derivative can be seen as a local differential operator which acts on coho¬ 
mology classes trivializing them [84], This fact has important consequences for symmetries and 
conservation laws associated with local variational problems generating global Euler-Lagrange 
expressions; see, e.g., [7, 25, 28]. 

As an immediate example of that feature, let (5(-) denote a cohomology class defined by 
a closed class in the variational sequence, and let e be a dynamical form locally variational: 
En+i{£) = 0 (i.e., locally, on an open set U,,, e = Ex^ = EniXi))- In the case q = n + 1, 
Theorem 4.10 implies that for every projectable vector field S on Y, Ljs^e = E„(J®SvJe), that 
is, the Euler-Lagrange form of the global Lagrangian A = which implies S{Ljs^E\J = 

6{ELjrsX,) = 0 even if the cohomology class of Ex^ is nontrivial. In particular, since E^jr^x, = 
En{Ljr^X^) = £'„( ExJ we see that Euler-Lagrange equations of the local problem defined 
by Aji-sAt are the same as the Euler-Lagrange equations of the global problem defined by 
A = J^EvlEx,- Summarizing, a problem defined by a family of local Lagrangians Ljr^X^ is 
variationally equivalent to a global one [84, 89]. 

Analogous results can be obtained at any higher degree in the variational sequence. Another 
case of particular importance, the degree n + 2, will be discussed below. 

5.5 Symmetries of Helmholtz forms and variational dynamical forms 

Invariance properties of classes in the variational sequence suggested to Krupka et al. [56, 57, 58] 
the idea that there should be a close relationship between the notions of local variational triviality 
of a differential form and invariance of its exterior derivative. Namely, ‘translating’ Theorem 4.3 
to the Takens representation, we get [57, 58]: 

Theorem 5.2. The Lie derivative of a source (n + k)-form e by a contact symmetry Z is locally 
variationally trivial if and only if the flow of Z leaves invariant the source {n + k + l)-form 
En+k (£■) • 

Applying this fact to (n -|- l)-forms we then get a rather surprising assertion: 

Corollary 5.3 ([57, 58]). Contact symmetries of Helmholtz forms transfer nonvariational dy¬ 
namical forms to locally variational ones. 

This result has an interesting application: namely, to a system of equations which is not 
variational as it stands, and even does not possess any variational multiplier one can find a va¬ 
riational system related by a contact transformation Examples showing this interesting feature 
are elaborated in the papers [57, 58]. 

Furthermore, one can obtain interesting results in study of global properties of variational 
problems. Let us consider a variationally trivial Helmholtz-like form p, i.e., let En-\-2{'>l) — 0) 
and denote the corresponding local system of dynamical forms by e^, i.e., let rj = En+iiefl- If 
a projectable vector field E on T is a symmetry of r] then, by the above, each of the local dyna¬ 
mical forms Si = Lja^s^ is locally variational (meaning that locally it is the Euler-Lagrange form 
of a Lagrangian). In fact, applying Theorem 4.10, case k = 2, using the Takens representation 
with rj = Rn+ 2 {[p\)-, and formula (4.6), we obtain 

0 — Lja^rj — Lf) — — 0 -|- 7?j2-i-2‘f^n+i('7 [p]) 

= En+lRn+l{J‘^'^v\ [p]) = En+lRn+l{[J^’Rv Iv]) = -E'n-|-lT( J®EyJ 7?), 

SO that En+i{Ljs^£^) = En+iI{J^^v\p) = 0, meaning that the system of local Euler-Lagrange 
forms Lja^e^ is variationally equivalent (in the sense that they have the same Helmholtz form) 
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with the global dynamical form e = Moreover, e is locally variational, and therefore 

locally it is the Euler-Lagrange form of a local Lagrangian [85]. 

We conclude that the variational problem defined by a system of local Euler-Lagrange 
forms St is variationally equivalent to the global variational problem defined by the Euler- 
Lagrange form e. Remarkably, in mechanics this ‘globalization procedure’ means to add to the 
local equations a variational force, ‘sewing the equations together’ on overlaps of their domains. 

5.6 Generalized symmetries generating Noether currents 

Representing the variational sequence by differential forms is of fundamental importance for 
understanding the geometry of Noether theorems. 

Let e = E\^ be a locally variational form on J^Y, and let a projectable vector field E on 
be a generalized symmetry, meaning that Ljs^e = 0. Recall that by the properties of the Lie 
derivative this means that En{Ljr^X^) = 0, hence locally Ljr^X^ = dnldu- On the other hand, 
by Theorem 4.10, Ljs’^s = E„(J^EyJe), hence = 0 , and from here we have locally 

= dn^i, where are local currents which are conserved on-shell (i.e., along critical 
sections). Denoting by e,, the usual canonical Noether current defined by Proposition 4.7, case 
q = n, i.e., e, = J^Ey\pdv\,+J'''^H\K, and putting= n,+e„ we can write J^Ev\e+dHieL-/3,) 
= 0. We call the (local) current — [3^ a Noether-Bessel-Hagen current. Thus Noether-Bessel- 
Hagen currents — /3t are local currents associated with a generalized symmetry (conserved 
along critical sections); in [28] it was proved that a Noether-Bessel-Hagen current is variationally 
equivalent to a global current if and only if 0 = 6n-i{^v\SniK)) £ iL^(W, ker d//). 

It is of interest to investigate under which conditions a Noether-Bessel-Hagen current is 
variationally equivalent to a Noether conserved current for a suitable Lagrangian. A Noether- 
Bessel-Hagen current ex, — /3t associated with a generalized symmetry of e = Ex^ is a Noether 
conserved current if it is of the form ex, with fi, a current satisfying L = 0 

(i.e., E is a symmetry of the Lagrangian At = At — dn^i for e) [87]. The reader can check that 
taking for an invariant Lagrangian instead of the canonical cnrrent J^E\0 a current, related to 
a Lepage equivalent different than the Cartan form 9, we get a current which provides the same 
conservation law, see [68]. We have the on shell conservation laws dniex, — Pl) = 0. It can be 
proved that ex, — Ljs'^gL, = ex,-dHiJ'i. closed but also exact on shell (for a detailed 

discussion see, e.g., [II]). In the following we shall give an example of application concerning 
symmetries of Jacobi type equations for a Chern-Simons 3D gauge theory. 

Example 5.4. Let us then consider the 3-dimensional Chern-Simons Lagrangian 

-^cs(^) = ^AidjAjji + —AiAjAm^ uiQ, 

where ojq is a 3-dimensional volume density, n a constant, while Ai = A^J^ are the coefficients 
of the connection I-form A = Aidx^ taking their values in any Lie algebra g with genera¬ 
tors J^. By fixing g = s[(2,]R) and choosing the generators Jy = whith Pauli matrices, 
we have [J^_l,Jv\ = p'^^ep^uJa and Tr(J^Jj,) = \ripv, with rj = diag(-I,I,I) and eoi 2 = I- 
Hence, we can explicitly write Acs(^) = ~ where = 

diA>;-djA>^ + eypAf Aj is the so-called field strength. Note that, while Chern-Simons equations 
of motion are covariant with respect to spacetime diffeomorphism and with respect to gauge 
transformations, the Chern-Simons Lagrangian, instead, is not gauge invariant. Indeed for 
a pure gauge tranformation generated by a vertical vector field ^ - acting on the gauge potential 
as where Vj is the covariant derivative with respect to a principal connection uj^ 

defined by A^ in a usual way through the introduction of new variables [I] - we have £^Acs(^) 
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Let + be a 0-cocycle of strong Noether currents for the Chern-Simons Lagrangian (here i 
refers to an open set in a good covering and in the following we shall use the identification 
'^cs(^) = -^cs(^) and use them indifferently); it is known [27] that if is a symmetry of the 
Chern-Simons dynamical form, the global conserved current 

ShJ^hAcSi + f'^v\PdvC’3.\cSi,-: 

is associated with the invariance of the Chern-Simons equations and it is variationally equivalent 
to the variation of the strong Noether currents + Sl- Generators of such a global current lie 
in the kernel of the second variational derivative and are symmetries of the variationally trivial 
Lagrangian 

Let now require to be a symmetry of the Euler-Lagrange class of the ‘deformed’ Lag¬ 
rangian Cjs^Xcs(L), but not of the Chern-Simons class En{Xcs{t')) (this means that is not 
a gauge-natural lift). Note that in this case >Cjs=Acs(^) needs not to be a local divergence, in 
general. By naturality and by Theorem 4.10 this means that we require 


Cjs^En{Ev\En{Xcs{>')) = 0, (5.2) 

i.e., j^E is a symmetry of a Jacobi type dynamical form (or Euler-Lagrange class of the deformed 
Lagrangian). Let now j^Ey be the generator of a gauge transformation. Since we know that 
Ej‘>SyXcs{i') = dn'J, it is clear that such a j^E is a Noether symmetry of the Lagrangian 

SyjE„(Acs(0 - duij - f'^v\Pdv\cs)'> 
this comes from equation (5.2) and from 


Cj»'^{Cjs’=Xcs{i^) - dH^Xcs) = - dH{^v\Pdv\cB))- 


By applying the above results concerning Noether-Bessel-Hagen currents to the ‘deformed’ Lag¬ 
rangian we can then state that along sections pulling back to zero the Jacobi type dynamical 
form (note that solutions of Chern-Simons equations are among such sections) we have that 

dH[j'^v\Pdv(C^^\cs-dH^v\PdyXc^)\ = 

i.e., the corresponding Noether-Bessel-Hagen current must be exact on shell which implies a con¬ 
straint for the connection 1-form A. 

5.7 Sharp obstruction to the existence of global solutions 
in 3D Chern-Simons theories 

We now explore the relationship between the existence of global critical sections and the existence 
of global Noether-Bessel-Hagen currents. We denote by [[-jldR a class in de Rham cohomology. 
As already recalled, the Noether-Bessel-Hagen current is variationally equivalent to a global 
current if and only if 0 = Sn-iC^vl^) ~ [[-\/Je]]dR G 77 ^r(^)- 

In the following we shall need to see certain cohomology classes in as the pull-back 

of cohomology classes in i.e., defined by closed differential forms on the base manifold. 

Let us then assume cr be a global section of tt : Y —>• W; for any global section we have of course 
a* OTT* = idj:/nj^(x)- Let furthermore H2^{Y) ~ 7r*(L7^j^(X)), i.e., the pull-back tt* : —)• 

H^^fY) is an isomorphism of cohomology groups. 

We need to compare the cohomology of the base manifold with the cohomology of the to¬ 
tal space, and while H^^{Y) ~ Hy^fY), in principle we do not have a similar isomorphism 
for the cohomology groups H2^{X). This is overcome by the peculiar structure of quotient 
space in the variational sequence: since for all sections = 0, we can correctly define 
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J’'cj*[[HyJe]]dR = [[J'’cT*(HyJe)]]dR; in other words, since 0^ consists of contact n-forms the pull¬ 
back a* factorizes through A^/0^. Furthermore, there always exists a closed form a G A^, e.g., 
a Lepage equivalent, which represents the cohomology class [[EyJe]]dR and which projects onto 
Hyje G A^/0^. Thus, if cr is a critical section, = 0 and the corresponding class vanishes 

in therefore, by our assumption on the nth cohomology groups [[EyJe]]dR = 0 and 

there exist global Noether-Bessel-Hagen conserved currents for all generalized symmetries S [28]. 

We can use the result above also conversely to define a cohomological obstruction to the 
existence of global solutions for a given problem. Let then cr be a section of Y over X. Let 
a be a closed differential form which represents [[SyJejjdR G the obstruction to the 

existence of global conserved quantities for conservation laws associated with the problem e and 
its symmetry S. Since a (as well as any other section) defines a projection a*: H2^{Y) i—)• 
if the class T^(T*[[a]]dR £ does not vanish, then there are no global solutions. 

In the case we drop the above condition on the isomorphism of cohomology groups, if the 
pull-back with a given section of the cohomology class [[SyJejjdR is not zero, no other homotopic 
section can be a critical section. Note that, in this case, the obstruction does depend on the 
section; J^(T*[[a]]dR may vanish for sections from one homotopy class, but not for those from 
another [28]. 

It is also important to note that this obstruction is not identically zero, in general. 

As an example, let us consider the classical 3-dimensional Chern-Simons theory. It is a classi¬ 
cal field theory for principal connections on an arbitrary principal bundle P over a 3-dimensional 
manifold M. It is well known that the equations can be written F = 0 (F the curvature); these 
equations admit solutions if and only if the cohomology class [[trT) 4 ]]dR (“first Chern class”) 
vanishes for an arbitrary (and, thus, for every) connection A. Thus, we will consider only princi¬ 
pal bundles P with group U(l) or C* (taking the trace of the curvature corresponds to “taking 
the determinant of P”). 

To apply our result, we have to formulate Chern-Simons theory in terms of the bundle of 
connections: the bundle of connections is defined as 

C := J^P/G^ P/G ^ X, 

the principal connections on P are in one to one correspondence with the (global) sections of the 
bundle C ^ X; Chern-Simons theory can then be formulated as a variational theory on J^C. In 
this case is it advantageous to lift the dynamical form of the problem to the de Rham complex: 
we can find a differential form S on J^C such that a section ct of C i—?■ X is critical if and only if 

JV(HJS) = 0, 

for all vector fields H on J^C. 

We have that the contact structure of J^P defines a connection on the principal bundle 
J^P I—)• C; its curvature is given by P = dAj/\dx^ (structure group U(l) or C*). This connection 
is universal in the sense that every principal connection on P i—)> X is “induced” by it via a section 
cr: X I—)• C, in particular Pg- = a*F. 

We can set now 

E := ttIqF a F, 

where S is, of course, the exterior differential of the Poincare-Cartan form. Notice that a pro- 
jectable vector field H is a symmetry of S if and only if (i(HjS) = 0. The corresponding class 
to [[HyJe]]dR for a generic e is now [[“:jS]]dR = 27r]'o((“:jP) AP). 

If the Euler-Lagrange equations admit global solutions, then P = dy and, thus, 

[[SjE]]dR = 27ri*,o((HjP) A dy) = 2d{7rl,{E\F A ^). 

If the equations do not admit global solutions then [[P]]dR 7^ 0. 
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Now, if Ai is a closed manifold, then there is a closed cohomologically nontrivial 1-form a 
such that a f\Frj and, hence, (7r*(a)) A T are cohomologically non trivial (Poincare duality). 

Since JV is parallelizable, we can easily find a vertical vector field such that ‘hj = 7r*(a) 
and 


[[$jS]]dR = 2<o(vr*(«)A.F). 

Thus, for closed 3-manifolds Chern-Simons theories admit solutions if and only if 

[[‘&JS]]dR = 0. 

For non compact manifolds this is, of course, no longer true; thus, for non compact 3-manifolds 
the question, whether the obstruction is sharp, depends on the existence of a suitable com- 
pactification, see also [88]. 
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